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> ■ Abstract 

. We study the big-bang nucleosynthesis (BBN) with the long-Hved exotic particle, 

called X. If the lifetime of X is longer than ~ 0.1 sec, its decay may cause non- 
thermal nuclear reactions during or after the BBN, altering the predictions of the 
standard BBN scenario. We pay particular attention to its hadronic decay modes 
5^ \ and calculate the primordial abundances of the light elements. Using the result, we 

derive constraints on the primordial abundance of X. 

Compared to the previous studies, we have improved the following points in our 
analysis: The JETSET 7.4 Monte Carlo event generator is used to calculate the 
spectrum of hadrons produced by the decay of X; The evolution of the hadronic 
shower is studied taking account of the details of the energy-loss processes of the nu- 
clei in the thermal bath; We have used the most recent observational constraints on 
the primordial abundances of the light elements; In order to estimate the uncertain- 
ties, we have performed the Monte Carlo simulation which includes the experimental 
errors of the cross sections and transfered energies. 

We will see that the non-thermal productions of D, "^He, ^He and ^Li provide 
stringent upper bounds on the primordial abundance of late-decaying particle, in 
particular when the hadronic branching ratio of X is sizable. We apply our results 
to the gravitino problem, and obtain upper bound on the reheating temperature 
after inflation. 



1 Introduction 



In modern cosmology, big-bang nucleosynthesis (BBN) is one of the most important sub- 
jects. In the standard scenario, neutrons freeze out from the thermal bath when the cosmic 
temperature is ~ 0.7 MeV and then the light elements (i.e., D, ^He, ^Li, and so on) are 
synthesized subsequently. As we will discuss in the next section, prediction of the standard 
BBN (SBBN) scenario is in a reasonable agreement with the observations. 

Predicted abundances of the light elements are, however, very sensitive to the cosmolog- 
ical scenarios. In particular, if we consider exotic cosmological scenarios based on physics 
beyond the standard model, theoretical predictions on the light-element abundances may 
be too much affected to be consistent with the observations. Thus, the BBN provides 
significant constraints on the new particles which change the cosmological evolution at the 
cosmic time t ~ 10~^ — 10^^ sec. If we consider physics beyond the standard model, there 
exist various candidates of such exotic particles. (Hereafter, we call such particle as X.) 

One example of the long-lived particles is the gravitino in supergravity theory PJ. 
Gravitino acquires mass from the effect of the supersymmetry breaking. In addition, 
its interactions are suppressed by inverse powers of the gravitational scale and hence its 
lifetime becomes very long (if it is unstable). In particular, for supersymmetric models 
with the gravitino mass ^3/2 ~ O{10^~^) GeV, lifetime of the gravitino becomes much 
longer than 1 sec and its decay may significantly affect the light-element abundances. 
(This is called the "gravitino problem.") Thus, the BBN provides substantial constraints 
on the properties of the gravitino and also on the cosmological scenarios. (For more details 
of the gravitino problem, see |2] and references therein.) 

In addition, even for the case where the gravitino is the lightest superparticle (LSP), 
the next-to-the lightest superparticle (NLSP) has long lifetime since the NLSP decays into 
its superpartner and the gravitino. The BBN imposes significant constraints on the case 
where the NLSP is the neutralino or the scalar-r Furthermore, moduli fields in the 

superstring theory are another candidates of X. Some of the moduli fields may acquire 
non-vanishing amplitude in the early universe. If so, their coherent oscillation may decay 
at a very late stage of the evolution of the universe. 

The exotic particles listed above are some of the famous examples and, if one considers 
particle-physics models beyond the standard model, there may exist long-lived particles 
which affect the BBN. Thus, it is important to study the BBN scenario with late-decaying 
particles. Such studies have been done by various groups IHl El El 13 UHl UH 1121 UHl UH 

unmniiniiiHj- 

In most of the previous studies (except [HI El |H1 1121) , however, hadronic decay modes of 
X were ignored although, for many of the candidates of the long-lived exotic particles, it is 
expected that the hadronic branching ratio is sizable. For example, even if X dominantly 
decays into photon (and something else), the hadronic branching ratio is expected to 
be non-vanishing since the (virtual) photon can be converted to the quark anti-quark 
pair. In this case, the hadronic branching ratio is estimated to be at least ~ ae^/Air ~ 
-|^g-(2-3) (^^j^]-^ Q,^^ being the fine structure constant), unless the hadronic decay mode 
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is kinematically suppressed. Of course, if X directly decays into quarks and/or gluons, 
hadronic branching ratio can be close to 1. If the massive particles decay into quarks 
or gluons during/after the BBN epoch, many mesons (mostly pions) and nucleons are 
produced. The emitted hadrons lose their energy via the electromagnetic interactions and 
scatter off the background nuclei. The emitted hadrons affect the BBN via two effects. 
One is the "inter-conversion" effect; if X decays at relatively early stage of the BBN (i.e., 

10^ sec), emitted hadrons may change the neutron-to-proton ratio. On the contrary, 
at the later stage of the BBN (i.e., t^lO^ sec), hadrodissociation processes are caused 
by energetic hadrons generated by the decay of X. Due to these effects, light-element 
abundances can be significantly affected. 

The purpose of this paper is to investigate the BBN scenario with the long-lived exotic 
particle using the currently available best knowledges on particle physics, nuclear physics, 
and astrophysics, paying particular attention to hadronic decay modes of X. We calculate 
the abundances of the light elements including relevant hadronic scattering processes (as 
well as photodissociation processes). Then, in order not to spoil the success of the BBN, 
we derive upper bounds on the primordial abundance of X. 

This paper is the full-length version of our recent letter JH] ■ This paper is organized 
as follows. In the next section, we briefly review the current status of the observations 
and SBBN. Then, in Section |H1 we give an overview of the decay of massive particles and 
its cosmological effects. In Section |3J we give a brief overview of the photodissociation 
process. We outline the hadronic decay scenarios in Sectional In Section IHl we introduce 
the formulations and computations of inter-conversion effects between background p and n 
by hadrons at earlier epochs in the hadron injection scenario. In Section [7| we discuss the 
destruction and production processes of light elements in hadrodissociation scenario. We 
also consider the non-thermal production processes of Lithium and Beryllium in hadronic 
decay scenario in Section |S| In Section El we compare the theoretical predictions with the 
observations in hadronic and radiative decay scenario for general massive particles. Our 
main results are shown in this section; if the reader is mostly interested in the resultant 
constraints, see this section (in particular, see Figs. EH — 1121) • Then, in Section ITUl we 
apply our results to the case of decaying gravitinos in supergravity. Section [TT] is devoted 
to the conclusions and discussion. 

2 Current Status: Observation and SBBN 
2.1 Current status of observations 

First we briefly summarize the current status of the observational light element abun- 
dances. The errors of the following observational values are at la level unless otherwise 
stated. 

The primordial value of the ratio wd/'t^h is measured in the high redshift QSO absorp- 
tion systems. (Here and hereafter, denotes the number density of the nucleus Ai.) 
Recently a new deuterium data was obtained from observation of the absorption system 
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at the redshift z = 2.525659 towards Q1243+3047 [23 including improved modeling of 
the continuum level, the Ly-a forest and the velocity structure of the absorption systems. 
The reported value of the deuterium abundance by using Keck-I HIRES, was relatively 
low, iriB/nnT^' = (2.42l[J:|) x 10"^ Combined with the previous data [21 123 EEl EI! , it 
is reported that the primordial abundance is given by *^ 

{n^/nnr^ = (2.78+°;^^) x 10-^ (2.1) 

which we adopt in this paper. (Here and hereafter, the superscript "obs" is used for the 
primordial values inferred by the observations.) We should make a comment on taking 
mean of the data. Although we think that this is reasonable as treatment of experi- 
mental data, it was pointed out that the five measurements have a large dispersion than 
expected 123- Since the purpose of our paper is to derive a conservative constraint, we 
have to care about possible systematic errors. Since non-thermal production of D leads 
to a severe constraint on the abundance of massive particles, the upper bound of D/H is 
important for us. The highest value among the five measurements is fl^ |2S1 12nj 

(nD/nH)°"^ = (3.98l°:f7) X 10-^ (2.2) 

Thus, when we derive the constraint on generic massive particles in Section |^ we will also 
show the result with adopting Eq. ()2.2p . 

The primordial abundance of *^He is inferred from the recombination lines from the 
low metallicity extragalactic HII regions. One obtains the primordial value of the ^He 
mass fraction Y by regressing to the zero metallicity 0/H — > for the observational data. 
Based on the reanalysis of Fields and Olive jHOl, which takes account of the effect of the 
Hel absorption, the primordial mass fraction is given by 

yobs^pQ) = 0.238 ± (0.002)stat ± (0.005)syst, (2.3) 

where the first and second errors are the statistical and systematic ones, respectively. On 
the contrary, Izotov and Thuan |3II reported a slightly higher value : 

F°^^(IT) = 0.242 ± (0.002)stat(±(0.005)syst), (2.4) 

where we have added the systematic errors following Refs. [S21 EHl Ell • Since there exists 
sizable difference between the results of two groups, we use two values F°^'^(FO) and 

yobsj^jrj-j_ #2 

*^Note that higher deuterium abundance in relatively low redshift absorption systems at z = 0.701 
was also reported: no/j^H = (2.0 ± 0.5) x 10"'* [SH]. Based on another independent observation of the 
clouds, however, it is claimed that the observed absorption is not due to D although there are still some 
uncertainties. Thus, we do not adopt the "High D" primordial abundance in this paper. 

#2Recently Olive and Skillman ^ re-analyzed the Izotov-Thuan data of 1998 EEl (and 2004 [SI) 
and obtained "^Hc abundance with larger uncertainties. They used only ''He lines to estimate the "^He 
abundance, electron density and temperature. This "self-consistent" approach has some merit but cannot 
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As for ''Li, it is widely believed that the primordial abundance of ^Li can be deter- 
mined using Pop II old halo stars with temperature higher than ~ 6000 K and with low 
metallicity. We use the most recent measurements by Ref. logio {iT-'^u/n-uY^^ = 

-9.66 ± (0.056)stat ± (0.06)sys, which corresponds to nru/riu = (2.19t|jj^) x 10-^°. It was 
claimed that there can be a significant dependence of nru on Fe abundance in the low 
metallicity region [SHI- In addition, assuming that this trend is due to the cosmic ray 
interactions, Ref. jlU] inferred that the primordial value is nri^i/riYi = (1.23lo;32) x 10^^°. 
This differs by a factor of two from the result given in ^T] . This suggests that the system- 
atic errors in both observations may be underestimated. We are afraid that more ''Li in 
the halo stars might have been supplemented (by production in cosmic-ray interactions) 
or depleted (in stars) [SHI- Since the precise determination of the primordial abundance 
from the observations is out of the scope of this paper, we conservatively adopt the center 
value given by [SZ] with larger uncertainties in this paper: 

logio [{nru/riuT''''] = -9.66 ± (0.056),tat ± (0.3)add, (2.5) 

which corresponds to nrj^^/nu = {2.19^1'^) x 10^^°. This is also justified from the point 
of the view of deriving conservative constraints; with the primordial abundance of ''Li 
given by [10] , discrepancy between the values of the baryon-to-photon ratio 1] determined 
by the SBBN and that by the observations of the cosmic microwave background (CMB) 
anisotropies becomes worse. (See the next subsection.) 

For ^Li, it is much more difficult to determine its primordial abundance since ^Li is 
much rarer than ''Li. Unfortunately, data is insufficient and ^Li abundance cannot be 
reliably determined. However, because it is generally believed that the evolution of ^Li is 
dominated by the production through the cosmic ray spallation (i.e., reactions of cosmic 
rays with the interstellar medium), we can set an upper bound on the ratio neu/nri^i. The 
models of the nucleosynthesis through the cosmic ray spallation were intrinsically required 
to simultaneously explain the whole observational Li-Be-B abundances |1T1021E31- On 
the other hand, recently it was claimed that the observational ^Li abundance in halo stars 
is too abundant from the point of view of the cosmic ray energy if ^Be is fitted by the 
model of the cosmic-ray metal |31]. Therefore, there seems to be some uncertainties in 
the models of the cosmic ray spallation. *^ In this situation, at least it would be safe to 

determine the temperature and density precisely as the authors admit. The standard approach adopted by 
Izotov and Thuan (and others) is to use the OIII Hues to estimate the temperature because it leads to much 
more precise determination. At present we think that it is premature to judge which approach is more 
adequate to estimate the ^He abundance. Therefore, we did not adopt it in this paper. For reference's 
sake, if we adopted the value given in Ref. |35| . the constraint on the abundance of the massive particle 
(or the reheating temperature after the inflation) would become milder by a factor of 4 ~ 5 compared to 
the result with Izotov and Thuan's value at short lifetime (^; 10^ sec). 

"^^^Recently, Suzuki and Inoue |45| pointed out other possibility of producing ®Li independently of the 
abundance of ^Be through a-a reactions induced by cosmic-ray a accelerated in structure formation 
shocks. However, it would be difficult to precisely predict the abundance of ^Li in the current version of 
their model. Therefore, it is premature to quantitatively discuss the abundance by using their model. 
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assume that ^Li abundance increases as the metalhcity increases. Today we observe only 
the ^Li to ^Li ratio in low-metalhcity ([Fe/H] < —2.0) halo stars 



(n6Li/n7Li)^^^° = 0.05 ± 0.02 {2a). (2.6) 

We take this value as an upper bound on the primordial value ofnei^-Jmu- (^SLi/'^'^Li)"'^'^ < 
(?T,6Li/r;,7Li)^'^^°. In our statistical analysis, we use the ratio neLi/nn in deriving the con- 
straints since, with the long-lived exotic particle X, this ratio can be calculated more 
reliably than neLi/riTLj. In particular, combining Eq. (j2.6p with Eq. (|2.5|) . we use the 
upper bound 

(mLi/nH)°"^ < (n6Li/nH)^"'° = (1.10l^;°°) x lO"" (2a). (2.7) 

For the constraint on ^He, we adopt the observational ^He abundance of the pre-solar 
measurements. In this paper, we do not rely upon any detailed models of galactic and 
stellar chemical evolution because there are large uncertainties in extrapolating back to 
the primordial abundance. According to such theories of the chemical evolution, the ^He 
abundance can decrease or increase after the BBN epoch. Therefore, "^He-to-H ratio" in 
itself can not be solely used for a constraint. Instead we adopt the present ratio of ^He to 
D, r3 2, as the upper bound on the primordial value. This is based on the following simple 
argument of the chemical evolution. Suppose that some astrophysical process destroys D 
and/or ^He as 

AnsHe = -Rsn-iue, (2.8) 
Ann = -i?2nD, (2.9) 

where AnsHe (A^d) is the change of the ^He (D) abundance and R3 {R2) is the probability 
of destruction of ^He (D). Then, the change of the ratio rs 2 = n3He/''^D is 

Ar3,2 = — r = ^ ^'^3,2- 2.10 

riD + AriD riD I - R2 

Since D is more easily destroyed than ^He, *^it is quite reasonable to assume 

R2 > R3, (2.11) 

which leads to Ar3^2 > 0. Thus, the ratio r3^2 is monotonically increasing function of the 
cosmic time. Note that, in order to derive this result, we only rely on the fact that D 

"^^Recently, using the Subaru Telescope, ^Li/^Li was measured in the metal poor subgiant HD 140283 
with the use of a high-S/N and high-resolution spectrum: (nen/rtH)^'''" < 5.0 x 10"^^ {2a) gT]- The 
subgiant HD 140283 is the metal-poorest among all of the objects which have been used to derive the 
bound on ^Li/''Li ([Fe/H] — —2.5). It is, however, premature and beyond the scope of this paper to judge 
its reliability. In addition, our purpose is to obtain a conservative constraints. Thus, we do not use this 
constraint. 

#5The binding energy of D is 2.2 MeV while the threshold energies of any destruction processes of ^He 
are larger. 
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is more fragile than ^He. Therefore the present ratio gives us an upper bound on the 
primordial value of r^^2- When we adopt the solar-system data jlHl, the ^He to D ratio is 
given by 

= 0.59 ±0.54 (2a). (2.12) 
We take this to be an upper bound on the primordial ^He to D ratio 

<2 < r%. (2.13) 
Naively it means the upper bound rf^2 < 1-13 (2cr). 



2.2 Current status of SBBN 

In this subsection, we briefly discuss the current status of the SBBN. In the recent years, 
there have been great progresses in the experiments of the low energy cross sections for 
86 charged-particle reactions by the NACRE collaboration [IH]- In the compilation, 22 
reactions are relevant to the primordial nucleosynthesis, and the old data were revised. In 
particular, 7 reactions of them are important for the most elementary processes generating 
nuclei with atomic number up to 7. Cyburt, Fields and Olive reanalyzed the NACRE 
data |5(Hl51j. They properly derived the la uncertainty and the normalization of the center 
value for each reaction. In addition, they also reanalyzed the four remaining reactions, 
using the existing data [S21 IHni IM] and the theoretical prediction (for one reaction) [SH] • 
Their efforts are quite useful for the study of the Monte Carlo simulation in BBN, and it 
was shown that their treatment is consistent with the other earlier studies adopting the 
results of NACRE jHElEZI. 

In our numerical study of the SBBN, we used the Kawano Code (Version 4.1) jHSj 
with some updates of the cross sections of the nuclear reactions. We use the center values 
and errors of the cross sections for 11 elementary nuclear reactions given in Refs. [301 
IHT] . (For the neutron lifetime, see Eq. ()5.6|) .) To systematically take into account the 
uncertainties, we perform the fltting including both the observational and theoretical 
errors which are obtained in Monte Carlo simulation. (See the Appendix in f2]-) In our 
analysis, we assume that the theoretical predictions of nD/nn, and log]^Q[(n7Li/nH)] 
obey the Gaussian probability distribution functions with the widths given by the la 
errors. Concerning the observational values, they are also assumed to obey the Gaussian 
probability distribution functions. Note that we consider two cases for Y"^^, i.e.. Fields 
and Olive (FO) given in Eq. (|2.3|) and Izotov and Thuan (IT) given in Eq. (j2.4j) . 

We calculated the abundances of the light elements as functions of baryon-to-photon 
ratio: 

ri - ^, (2.14) 

where ub and are number densities of the baryon and photon, respectively. The results 
are plotted in Fig.^ As one can see, theoretical predictions become more or less consistent 
with the observational constraints when 77 ~ 6 x 10^^°. 
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Figure 1: Abundances of the light elements as functions of rj. The solid lines are the center 
value while the dotted lines show the theoretical uncertainties. Observational constraints 
are also shown. 
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Figure 2: variable as a function of 1] for SBBN with three degrees of freedom. For 
the constraint on Y, we used Fields and Olive's result (solid) and Izotov and Thuan's 
(dashed), which are given in Eqs. (j2.3|) and (j2.4j) . respectively. The shaded band (vertical 
solid lines) shows the baryon-to-photon ratio suggested by the WMAP at the la (2cr) 
level. 



For a more precise determination of 1], we calculate the x^-variable as a function of 1], 
and the result is shown in Fig. |21 The solid line (dashed line) is for the case of Fields and 
Olive (Izotov and Thuan). From this figure, we see that the prediction of the SBBN agrees 
with the observation of ^He, D, and ''Li at 95 % CL. In addition, we obtain the baryon 
to photon ratio at 95 % CL. as t]^^^^^) = 5.85;^]^ x 10-^° {^.QOtlH x 10-^°) using the 
value of Y in Fields and Olive (Izotov and Thuan). Since the baryon to photon ratio is 
related to the baryon density parameter as Qsh"^ = 3.67 x lO^r/, we obtain, at 95 % C.L., 

2 _ / OmU+lfofi (Fields and Olive) 
^^""^ ~\0.021410S (Izotov and Thuan) ' ^^'^^^ 

where h is Hubble parameter in units of 100 km/sec/Mpc. 

We also plot the value of r] reported by the WMAP collaborations in observations of 
the cosmic microwave background (CMB) anisotropies 123; which is approximately given 
by 

r/ = (6.1±0.3) X 10"^°, (2.16) 
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where we have adopted the shghtly larger error for the lower bound. The shadowed band 
in Fig. 121 represents the baryon to photon ratio at Icr. The vertical solid lines are their 
constraints at 2a . From this figure, we find that SBBN is consistent with the CMB 
observation. 

Under these circumstances, comparing the predictions of the BEN computations with 
observations, we can constrain the non-standard scenario such as the radiative decay or 
the hadronic decay of long-lived massive particles. 

Here we should mention that the consistency between theoretical predictions and ob- 
served abundances in SBBN or between CMB and SBBN is partly because we have adopted 
rather large systematic errors for the observed abundances of ^He and ^Li. In fact, if we 
adopted smaller systematic errors reported in the original papers, we would be confronted 
with difficulty that the 77 inferred from D and CMB disagree with that from ^He and 
''Li (e.g. see, Ref. [SHI)- -^^^ we assume large systematic errors and that SBBN is 
consistent because the purpose of the present paper is to derive conservative constraints 
on the massive particles with hadronic decay mode. 

3 Overview 

Before going into the detailed discussion of the BBN with long-lived particle X, we give an 
overview of the cosmo logical scenario we consider, and define parameters which are used 
in our analysis. 

3.1 Production 

In this paper, we consider a scenario where a massive particle X, with mass mx and decay 
rate Fx, has non- vanishing number density at the early universe. First, we consider the 
production of X in the early universe. Throughout this paper, we consider a situation 
where the X particle is somehow produced in the early universe. Production mechanism 
depends on the property of X. For example, if X is a particle, like gravitino, it can be 
produced by scattering processes of the thermal particles. In addition, it may be also 
produced by the decay of other particles. Moreover, condensation of some (exotic) scalar 
field may play the role of X. In such a case, non-vanishing initial amplitude of the scalar 
field provides non-vanishing number density of X at the late stage of the evolution of the 
universe. 

In order to perform our analysis as model-independent as possible, we do not specify 
the production mechanism of X. Indeed, constraints we will obtain depends only on the 
relic density of X (before it decays). In order to parameterize the number density of X in 
the early universe, we define the "yield variable" 
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which is defined at the time t -C . Here, nx is the number density of X while s is the 
total entropy density of the universe. Notice that, as far as we can neglect the entropy 
production, Yx is a constant when t -C ^x^. 

3.2 Decay 

In this subsection we discuss the decay of massive particles and its cosmological effects. 
The overview is schematically presented in Fig. 01 

3.2.1 Decay processes 

In studying the effects of X on the BBN, we classify the the decay process of X into two 
categories; radiative and hadronic decays. These decay processes cause different types of 
reactions and it is necessary to take account of both processes. Branching ratios decaying 
into radiative and hadronic particles depend on the model. In order to perform our analysis 
as model-independent as possible, we define the "hadronic branching ratio" 

jj X— »hadrons /o (-)\ 

= , 

J- X 

where Fx is the decay rate of X and Fx^hadrons is the hadronic decay width of X. 

If X may directly decay into colored particles, Bh may become close to 1. In addition, 
even when X primarily decays into photon (and other non-hadronic particles), is 
expected to be non-vanishing since the quark-anti-quark pair can be attached at the end 
of the (virtual) photon line. 

For example, for the case where the unstable gravitino t/'^ plays the role of X, which is 
one of the most well-motivated cases, the gravitino may directly decay into gluon-gluino 
and/or quark-squark pairs. (See Fig. HI) If the decay rate of these modes are sizable, Bh 
becomes close to 1. If these hadronic decay modes are kinematically blocked, however, the 
gravitino may primarily decay into photon and the neutralino x^: ^ 7 + If this is 
the only kinematically allowed two-body decay process of the gravitino, Bh becomes much 
smaller than 1. However, even in this case, it is also expected that Bh is non- vanishing 
since the decay mode ipfi ~^ Q + Q + (with q and q being the quark and anti-quark, 
respectively) exists. (See Fig. 0) In this case, Bh is expected to be Bh ~ lO"*^^"^-', since 
the process is three body process so the branching ratio is suppressed by the factor aem/4vr. 
To be more quantitative, we considered the case where the lightest neutralino is purely the 
photino 7, the superpartner of the photon, as an example. We calculated the branching 
ratio for the process ipfj. q + q + 'y (where q here corresponds to u, d, s, c, and b). In our 
calculation, the Feynman amplitude is calculated by using the HELAS package then 
the phase-space integration for the final-state particles is done with BASES package [UU] . 
Results for several cases are shown in Tabled As one can see, Bh in this case is indeed 
]^g-(2--3)^ (We have checked that Bh is insensitive to the gravitino mass as far as the mass 
difference ^3/2 — fri;^/ is fixed.) 
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Figure 3: Flow-chart of the hadronic decay of massive particles. 
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Figure 4: Feynman diagrams for the decay processes if)^ ^ g + g and ij)^ Q + Q, where ip^, 
g, g, q, and q are the gravitino, gluon, gluino, quark, and squark, respectively. Here, the 
black blob represents the vertex originating from the gravitino-supercurrent interaction. 




Figure 5: Same as Fig. but for the radiative decay modes. (Here, 7 is the photon while 
X*^ is the neutralino.) 



12 



^3/2 ^3/2 - ^7 



100 GeV 
100 GeV 
100 GeV 
300 GeV 
300 GeV 
300 GeV 



20 GeV 
40 GeV 
60 GeV 
20 GeV 
40 GeV 
60 GeV 



0.004 
0.008 
0.010 
0.003 
0.008 
0.010 



Table 1: Hadronic branching ratio of the gravitino for several values of the gravitino mass 
m3/2 and the mass difference between the gravitino and the photino — Here, the 
lightest neutralino is assumed to be the pure photino and all the superparticles other than 
the photino are assumed to be heavier than the gravitino. The minimal possible value of 
the invariant mass of the quark anti-quark system is set to be 10 GeV. 



We also introduce several parameters in order to characterize the decay of X. First, we 
assume that each primary parton jet has the energy E-^ct- (For example, when X decays 
into q-q pair, Ej^t = ^fnx-) This parameter is used when we study the hadronization 
processes with JETSET event generator. On the contrary, for the decay process of X with 
energetic photon in the final state, we define E!:^^ which is the energy of the emitted photon. 
In addition, in some case, invisible particle may be also emitted. For example, when the 
gravitino plays the role of X, some fraction of the energy is carried away by the LSP which 
we assume is the neutralino. Thus, we define E'vis, which is the (averaged) energy emitted 
in the form of the "visible" particles. As we will discuss, number of high-energy photons 
produced in the electromagnetic cascade process is proportional to i?vis- 

3.2.2 Radiative decay 

In the radiative decay the massive particles decay into photons (and/or electrons, and 
so on). The energy of the emitted particles can be as large as the mass of the parent 
particles; with the high energy primary particles, electromagnetic showers are induced 
and energetic photons are recursively produced in the shower (see the right branch in 
Fig. 01). Some of soft photons produced secondarily in the shower induce destruction and 
production processes of the light elements (D, ^He, '^He, ^Li and ^Li). For the details, see, 
for example, Ref. jnHHCH. 

In the case of radiative decay, the most important process that determines the destruc- 
tion rate of the light elements is "photon-photon process" where the high energy photons 
scatter off the background photons into electrons and positrons. Since the number of the 
background photons is about 10^° times larger than that of electrons, the photon-photon 
process, if it occurs, quickly thermalizes the high energy photons and the destruction of 
light elements does not take place frequently. However, since there exists the threshold of 
the photon energy [Eth ~ me/22T), soft photons with energy less than the threshold can 
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destroy the light elements. Comparing the binding energies of the light elements with E'th, 
we can see that the photodissociation processes of D and ^Hc may become effective when 
the temperature becomes lower than T^O.Ol MeV and 0.001 MeV, respectively. There- 
fore, roughly speaking, the constraint on the radiatively decaying particle comes from the 
destruction of D for hfetime less than 10^ sec and from the overproduction of D and ^He 
due to the destruction of ^He for hfetime longer than 10^ sec. 

3.2.3 Hadronic decay 

When quarks or gluons are emitted in the decay of the massive particles, they firstly 
fragment into a lot of hadrons and form hadronic jets. As a result, many high energy 
mesons and nucleons are injected into the cosmic plasma. 

At earlier epochs (t ^ 100 sec) the high energy mesons and nucleons lose their energy 
very quickly through electromagnetic interaction. Thus, they are completely stopped and 
reach to the kinetic equilibrium. Thus, the emitted hadrons do not directly destroy the 
light elements. After energy loss, they scatter off the background p or n through the strong 
interaction with their threshold cross sections. Then, they inter-convert the background 
p and n each other, even after the normal freeze-out time of the n/p ratio of the weak 
interaction. Since the typical mean free time for strong interaction is C(10"^) sec, only 
mesons with relatively long lifetimes such as tt^ and K^'"^ and nucleons (p, n, p, n) can 
cause the p-n inter-conversion. Since, at T ^ 1 MeV, the proton is more abundant than the 
neutron, the conversion from p to n takes place more frequently than its inverse process, 
and hence the hadron injection extraordinarily tends to increase the ratio n/p. As a result, 
the produced ^He and D would increase in the hadron injection scenario compared to the 
SBBN case. 

At later epoch (t^ 100 sec), mesons decay before they interact the background nucle- 
ons, and hence they become cosmologically irrelevant. On the other hand the emitted high 
energy protons and neutrons can scatter off the background p, n and "^He (which is synthe- 
sized in BBN). Since the energy loss due to the electromagnetic interaction is insufficient, 
the high energy p and n interact with background hadrons before they lose energy, and 
produce secondary hadrons through elastic and inelastic collisions. Such hadronic inter- 
actions occur successively and evolve hadronic showers. During evolution of the hadronic 
shower, a lot of "^He's are destroyed by the inelastic collisions, and D, T and "^He are 
produced from the ^He dissociation. Then the energetic T and ^He scatter off the back- 
ground ^He, and produce ^Li and ^Li. Since ^He is much more abundant than the other 
light elements, those non-thermal production of D, ^He, ^Li and ^Li drastically changes 
the prediction of SBBN. Thus, the severe constraint is imposed on the hadronic decay. 
Contrary to the radiative decay, non-thermal production due to ^He dissociation is even 
important at 10^ sec when the high energy photons quickly lose their energy by the 
photon-photon process and cannot destroy ^He. 

Here, we remark that the almost all of the energy of the primary hadrons are transferred 
to the electrons (positrons) and photons through the electromagnetic energy loss processes. 
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and decays of mesons and heavy charged leptons. Then the energetic photons and electrons 
cause electromagnetic showers. In this sense, the hadronic decay also has the same effect 
as the radiative decay. This is indicated by dashed arrows in Fig. El 

In order to take account of the non-standard processes, we have modified the Kawano 
Code (Version 4.1, with the updates of the nuclear cross sections). In particular, we have 
developed new subroutines which deal with the photodissociation and hadrodissociation 
processes and implemented them into the Kawano Code. Details will be explained in the 
following sections. 

3.3 Comparison with previous works 

Before going into the details of our analysis, it would be relevant to compare our analysis 
with previous ones. As we mentioned, once the hadronic decay occurs, the BBN processes 
are affected by the inter-conversion and hadrodissociation processes. 

The effects of the inter-conversion was first studied in whose results were also 
applied to some of the topics in the early universe in l2Uj . In these works, however, 
it is simply assumed that all of emitted hadrons are effectively stopped even when the 
cosmic time is longer than t ~ 10^ sec. Because the energy loss process becomes inefficient 
at t > 10^ sec for proton and t ^ 10^ sec for neutron, as we will show in detail in Sectional 
their assumption becomes inappropriate at low temperature. In our analysis, we carefully 
reconsider study the stopping of the hadrons through the electromagnetic interactions 
with background plasma. (Such stopping processes are also important for the study of the 
hadrodissociations.) On the other hand, effects of the hadrodissociation processes were 
studied in [8 . This study is, however, based on theoretical and experimental information 
which can be improved with our current knowledge. Thus, in our study of the evolution 
of the hadronic showers, the basic framework is the same as that used in [8^ , but there are 
several important modifications (see the following sections). 

After these pioneering works, there have been various theoretical, experimental and 
observational progresses to study the BBN scenario with long-lived exotic particles. First, 
with the progresses in the high-energy experiments, now we have better information on 
the hadron fragmentation processes. In particular, we use the JETSET 7.4 Monte Carlo 
event generator [2Tj to estimate the distributions of the nucleons and mesons produced 
by the hadronic decay of X. In addition, we have more experimental data of the hadron- 
nucleon cross sections and energy distributions of the hadronic particles generated by the 
hadrodissociation processes. With these improvements, we can perform a better study 
of the evolution of the hadronic shower. Furthermore, observational constraints on the 
primordial abundances of the light elements are also improved. 

In summary, the most important improvements which are made in this paper are as 
follows, (i) We carefully take into account the energy loss processes for high-energy nuclei 
through the scattering with background photons or electrons. In particular, dependence on 
the cosmic temperature, the initial energies of nuclei, and the background ^He abundance 
are considered, (ii) We use available data of cross sections and transfered energies of 
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elastic and inelastic hadron-hadron scattering processes as much as possible, (iii) The 
time evolution of the energy distribution functions of high-energy nuclei are computed 
with proper energy resolution, (iv) The JETSET 7.4 Monte Carlo event generator j2l] 
is used to obtain the initial spectrum of hadrons produced by the decay of X. (v) The 
most resent data of observational light element abundances are adopted, (vi) We estimate 
uncertainties with Monte Carlo simulation which includes the experimental errors of the 
cross sections and transfered energies, and uncertainty of the baryon to photon ratio [22] ■ 

4 Photodissociation 

We are at the position to discuss various processes induced by the hadronic (as well as 
the radiative) decay of the late-decaying particle X. In calculating of the abundances 
of the light elements, we take account of two types of dissociation processes of the light 
elements; one is the photodissociations induced by the energetic photons and the other is 
the hadrodissociations by the hadrons. Importantly, even if we consider hadronic decay 
modes of X, kinetic energy of the hadrons are eventually converted to radiation via the 
scattering processes. Thus, even in the case of the hadronic decay mode, it is important 
to consider the photodissociation processes. In this section, we first discuss the simple 
photodissociation reactions of the light elements. 

Once the late-decaying particle X decays in the thermal bath, most of the (visible) 
energy released by the decay of X is eventually converted to the form of the photon for the 
situation we are interested in. Then, the electromagnetic cascade processes are induced. 
In order to study the photodissociation processes of the light elements, it is necessary to 
understand the spectrum of the high energy photon generated by the cascade process. 
In our study, we have calculated the photon spectrum taking account of effects of the 
following processes: 

• Injection of the high energy photon from the radiative decay of X 

• Double photon pair creation (7 + 7bg — e+ + e~) 

• Photon-photon scattering (7 + 7bg 1 + i) 

• Compton scattering off thermal electron (7 + CgQ — > 7 + e~) 

• Inverse Compton scattering off background photon (e^ + 7bg ^ + 7) 

• Pair creation in background proton (and aec) (7 + Pbg ^ + e~ + p) 

For details of the calculation of the photon spectrum, see Appendix El 

One important point is that the energy distribution of the photon in the electromag- 
netic shower is mostly determined by the total amount of the injected energy and is 
insensitive to the detail of the primary spectrum of the injected high-energy particles. 
Thus, the photon spectrum depends on the temperature, number density of X, decay rate 
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of X, and the total visible energy released by the single decay of X which we call E^is', 
once these parameters are fixed, the photon spectrum is determined. 

In addition, if the temperature is high enough, pairs of charged particles like /x^/i^, 
TT+TT^, and so on, may be produced by the photon-photon scattering. Such pair-production 
processes, however, do not significantly change the photodissociation rates since photon 
spectrum for the photon energy relevant for the photodissociation processes is determined 
by the Compton scattering and the e+e~ pair creation in the nuclei. Thus, for the study 
of the photodissociation, we can neglect the pair production of charged particles heavier 
than the electron. If the hadronic particles are pair-produced, however, it may provide 
new sources of the hadronic particles. Such effects will be considered in the next section. 

Effects of the photodissociation is taken into account by including the following terms 
into the Boltzmann equations describing the evolutions of the light elements: 



dt 



- photodiss 



J2 / (th) dE^aA,^A,{Ej)f^{E^) 
+ J2^A, / dE^aA,-.AXEj)fj{E^), 



(4.1) 



where GAi^Aj is the cross section for the process A, + 7 — + ■ ■ ■ with the threshold 
energy -E^**^-*, and the summations are over all the possible reactions. (In this paper, Ai is 
used for light elements, i.e., ra, p, D, T, ^He, ^He, ^Li, ''Li, and ^Be.) The above equation 
for z-th nucleus can be also written as 



dnA. 
~dt 



where 



Ep (photodiss) I Y^r 



(photodiss) 



riA, 



photodiss 



T (photodiss) 
Aj ^Ai 



= L(t^)'^El^A,^Aj{E^)f-f{E^)- 



(4.2) 



(4.3) 



Thus, the BEN reactions with the radiatively decaying particle can be followed once the 

J--J-- m (photodiss) 1 

quantities 1 Xj^Ai known. 

In our analysis, we have calculated r^^*^^'***^'* for all the relevant processes. In order 
to derive the constraint on the primordial abundance of X, we calculate the abundances 
of the light elements up to ^Li. All the processes included in our analysis are listed 
in Table |21 For the photodissociation of D, we use the cross section in the analytic 
form. For the other processes, the cross sections are taken from experimental data. For 
the photodissociation cross sections of ^Be, we could not find experimental data. Thus, 
we use the photodissociation rates of ^Li for corresponding dissociation processes of ^Be; 



T (photodiss) 
7Be 



T (photodiss) 
7Li^yl 



We have checked that, in deriving the constraints on the properties 
of X, ^Be does not play a significant role. Thus, our results are not significantly affected 
by this assumption. 
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Reaction 



Error Reference 



7 + D ^ n + p 
7 + T ^n + D 
'f + T—*p + n + n 
7 + ^He ^ p + D 
7 + ^He — > p + p + 
7 + ^He ^ p + T 
7 + ^He n + ^He 
7 + ^He — i> p + n + D 
7 + ^Li — s> anything 
7 + ^Li -> n + ^Li 
7 + ''Li — > anything 
7 + ^Be ^ p + ^Li 
7 + ''Be anything 



6 % inn 
14% inainsi 



7% insi 

10% lEl 

15% IMl 

4% inu 



5% inniEii 



14% ins! 

4% IHH! 

4% IMl 

9% HOI 



Table 2: Photodissociation reactions included in our analysis. We also write down the 
errors of each cross sections which are used in our Monte Carlo analysis. 



Since the abundance of ^He is large, the photodissociation of ^He may significantly 
change the abundances of nuclei with atomic number A < 4. Thus, in considering the 
photodissociation of nuclei with A < 4, we specify the final-state in order for precise 
calculation of the abundances of the light elements. The dissociation of Li, on the contrary, 
do not change the abundances of the nuclei lighter than Li because of the smallness of the 
abundances of ^Li and ''Li. Thus, for most of the photodissociation processes of Li, we do 
not specify the final-state particles. The only exception is the process ^Li + 7 — n + ^Li; 
this process may be important for the calculation of the ^Li abundance; In our analysis, 
the effect of the ^Li production through this process is properly taken into account. 

So far, we have discussed the scatterings of the photons with the background nuclei. 
Importantly, since the photons are energetic, the final-state particles produced by the 
photodissociation processes may acquire sizable energy and participate in other class of 
non-thermal production processes of the light elements. In particular, energetic T and ^He 
may scatter off the background a to produce ®Li. Since the observational upper bound 
on the primordial abundance of ^Li is very stringent, non-thermal production of ^Li may 
impose significant constraint on the properties of X |T^ ll4j. Such non-thermal production 
of ^Li will be discussed in Section jS] in detail. 

5 Hadronic Decay of Massive Particle 

Now, we consider hadronic decay of X and its effects on the light-element abundances. 
For this purpose, it is important to understand how the partons emitted by the decay of X 
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Vs (GeV) 



Figure 6: Charged-particle multiplicity {Nch) per two hadronic jets as a function of the 
center of mass energy ^/s [7T]. The solid line is the case with e~^e~ annihilation, while the 
dashed one is for the case with p + p{p) collisions. For the p + p{p) case, we use PYTHIA 
Monte Carlo event generator and do not include the single diffraction. 

are hadronized and how the hadrons propagate in the thermal bath. Thus, let us discuss 
these subjects in this section. 

5.1 Hadronization 

We first consider the hadronization processes. Since we mostly consider the cases where 
the mass of X is larger than the QCD scale, X primarily decays into the quarks and/or 
gluons for the hadronic decay process. For the cosmic temperature we consider, however, 
the time scale for the hadronization is much shorter than that for the cosmic expansion. 
Thus, the partons emitted from the decay of X are instantly hadronized and fragment into 
the mesons and nuclei (vr^, 7r°, K^, g, n, p, , and so on). In studying the evolution 
of the cascade processes, which will be discussed in the next section, those energetic nuclei 
and mesons (in particular, p and n) become the primary particles which cause the hadronic 
shower. 

In our study, the hadronization processes are followed by the JETSET 7.4 Monte 
Carlo event generator [21], which computes the fragmentation of the hadrons from the 
partons. For the decay process of X in which colored particles are emitted, we used the 
JETSET event generator to estimate the spectrum of the proton, neutron, pion, and so on. 
Importantly, predictions of the JETSET package well agree with the experimental data. 
Thus, in our analysis, we expect that the uncertainties arising from the hadronization 
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Figure 7: Spectra of the mesons (7r+ + vr^, + , and K^) produced by the e^e~ 
annihilation process as functions of their kinetic energy -Ekm- The center-of-mass energy 
is taken to be ^/i = 91.2 GeV. 

processes are small enough to be neglected. 

To demonstrate the agreements of the results from the JETSET package with exper- 
iments, in Fig. ini we plot the averaged charged-particle multiplicity (A'ch) (which is the 
total number of the charged hadrons) as a function of the center-of-mass energy of the 
qq system -^i. From Fig. IHl we find that the predicted {N^h) well agrees with the 
experimental values. In addition, as we can see from the figure, (A^^ch) increases as y/s 
increases. 

In addition, in Figs. [Z|and|Hl we plot the spectra of the mesons (tt^ + vr", + K~, 
and KD and baryons (n + n, and p+p), respectively, as functions of their kinetic energy in 
the center-of-mass frame E^^^^K Here we take the total energy of the jets to be ^/s = 91.2 
GeV, for which we see that the typical energy of the produced hadrons are (9(10) GeV. 
In studying the effects of the hadronic decay of X, we calculate the spectra of the primary 
hadrons with the relevant total energy of the jets which is determined by the mass of X 
and the decay mode. For mx ~ O(IOO) GeV — (9(100) TeV, we have found that the 
typical energy of the produced hadrons is about 0{1 — 100) GeV 

In Fig. ini we also plot the averaged number of the produced hadron per two hadronic 
jets as a function of the center-of-mass energy y/s. From this figure, we can see that the 
hadrons directly emitted from the decay of X is mostly pions. In addition, in Fig. ITUl 

"^^Here (Nch) is defined as the value after the decay of the Ks and A° since their lifetimes (r^o = 
0.89 X lO"""^" sec and tao = 2.63 x 10~^° sec, respectively |7T]) are much shorter than the cosmological 
time scale we are interested in. For details, see the subsection 15. 21 
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Figure 8: Spectra of the baryons ((a) n+n and (b) p+p) produced by the e~^e annihilation 
process as functions of their kinetic energy -Ekin- The center-of-mass energy is taken to be 
v^ = 91.2 GeV. 
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Figure 9: Number of the hadrons produced by the e^e~ annihilation process as functions 
of the center-of-mass energy ^/s. The dotted, the short dashed, the thin sohd, the dot- 
dashed, and the long dashed lines are tt^ + tt", + K~ , K^, P+P, and n + n, respectively. 
We also plot the charged-particle multiplicity Nch by the thick solid line. 




Figure 10: Number of the hadrons produced by the p + p{p) scattering process. Meanings 
of the lines are the same as Fig. IHl 
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we plot the averaged numbers of the produced hadrons in p + p{p) colhsion. We see 
that the number of the secondarily-produced nucleon-anti-nucleon pairs is small for the 
center-of-mass relevant for our study v^<C(10) GeV. Therefore, we ignore them in this 
paper. 

In considering the hadronic processes in our analysis, we use the hadron fragmentation 
obtained for the process e+ + ^ q + q with ^/s = 2Ejet as a primary spectrum of the 
hadrons generated from X. Here, Ej^t is the energy of the primary jet, and will be related 
to rrix later. We have checked that our results are insensitive to the Lorentz and color 
structure of the qq system as far as the value of ^/s is fixed. 

As we mentioned in the previous section, there is another possible production process of 
hadrons via the scattering of the high energy photons emitted from X with the background 
photons. In particular, the lightest charged mesons, vr^, can be generated from the process 
7 + 7BG TT^ + TT". Such mesous may contribute to the inter-conversion processes between 
the proton and neutron, which will be discussed in Section IHl 

For the high energy photon with energy E-y, the center-of-mass energy of such process 



is y/s ~ sJ?)E^fT. Thus, pair-creation rate of the pion becomes sizable only when the 

temperature is high enough. Importantly, if the pair creation of the charged pions are 
effective, electron-positron (and other lepton) pair can be also produced without kinemat- 
ical suppression. Thus, we estimate the number of the pions produced by the radiative 
decay of a single X as 



where (T^_|_^^^+_|_^- and a^^^^i+^i- are cross sections of the pion and charged lepton pair 
creation processes, respectively. Those cross sections are evaluated at s = 2>E^^T. For the 
explicit formulae for these cross sections, see Here, we only take account of the 

charged pion production by the primary photons directly emitted from the decay of X. 
Since the energy of the photon is drastically reduced after several steps of pair creations 
(and inverse Compton scattering of electron and positron), the pair-creation processes are 
mostly induced by the photons directly emitted from the decay process of X. In addition, 
if the energy of the primary photon become large, hadrons heavier than the pions may 
be produced by the photon-photon scattering processes. However, heavier hadrons mostly 
decay before inducing the inter-conversion processes. (Exceptions are kaons for some cases. 
Since the effects of the kaons are not included in discussing the inter-conversion effects, we 
do not consider the pair creation of the charged kaons here. For details, see Section In 
addition, productions of heavier hadrons are kinematically suppressed. Thus, we neglect 
the productions of those heavier hadrons. 

5.2 Various time scales 

Once energetic hadrons are emitted into the thermal bath, which consists of 7, e^, and 
nucleons, hadrons scatter off various background particles and induce cascade processes. 




^7r± 




(5.1) 
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Thus, in order to study the propagation of the hadronic particles in the expanding universe, 
it is important to understand the time scales of various processes; in particular, there are 
three important processes: (i) hadronic scatterings, (ii) decay (for unstable particles), and 
(iii) scatterings with thermal plasma through electromagnetic interactions. For the cosmic 
temperature we consider, the expansion rate is much smaller than the rates for the above 
processes and hence we can neglect the cosmic expansion in studying the effects of the 
processes (i) — (iii). Time scales of the electromagnetic processes will be discussed in the 
next subsection, and here we first consider (i) and (ii). 

First, we discuss the typical time scale of hadronic scattering processes. The interaction 
rate for the scattering processes between the projectile hadron Hi and the background 
nucleon Aj through Hi + Aj ^ + ■ ■ ■ is estimated as 



:.(4.4xlO-sec)-VA.f^) f^Spi^) (J_\\ (5.2) 
V ^ ^ V 6 / I 40mb ) VMeV/ ' ^ ' 

where, in this paper. Hi is used for hadrons. Here, is the velocity of iJj, ua is the 
number density of the nucleon species Aj, rjiQ = rjx 10^°, and /a = nA /ub with ub being 
the total baryon number density. For nucleon-nucleon collision processes, experimental 
data suggest that the product of (3 and a^j^^^ is approximately constant (~ 40 mb) for 
large range of the beam energies (see Figs. ITTl and IT^ . Thus, in the following discussion, 
we sometimes use 40 mb as a typical value of the cross section a^".^^^, although we use 
energy-dependent experimentally measured cross sections for our numerical calculations. 

The rate for the hadronic scattering process should be compared with the decay rate 
and also with the stopping rate in the thermal plasma. Among these two, we first consider 
the decay of the hadrons. Eq. (|5.2j) shows that the typical time scale for the hadronic 
scattering processes is longer than O{10~^) sec for T ^ 1 MeV. Thus, hadrons with lifetime 
longer than ~ 10^*^ sec participate in the hadron-hadron collision processes after the BBN 
starts. Thus, hereafter, we only consider such relatively long-lived mesons (tt^, K^, and 
Kl) and baryons (p, p, n, and n), whose lifetimes are |7Tj 

r^± = (2.6033 ± 0.0005) x 10^* sec, (5.3) 

tk± = (1.2384 ±0.0024) x 10~^sec, (5.4) 

Tj^o = (5.17 ± 0.04) X 10"^ sec, (5.5) 

r„ = 885.7 ±0.8 sec. (5.6) 

Other hadrons have very short lifetimes and decay before scattering off the background 
nuclei. 

#'^For example, lifetimes of 7r°, K^, and A° are r^o = 8.4 x 10^^'^ sec, r^o = 0.89 x 10^^", and tao — 
2.63 X 10^^° sec, respectively [7T]. 
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Since the lifetime of the neutron is relatively long, it is necessary to see if the non- 
thermally produced neutrons may decay before causing the hadrodissociation processes. 
As will be discussed in the next subsection, at low temperature T ^ me, energy-loss 
processes via the electromagnetic interaction becomes ineffective for non-relativistic neu- 
trons. Then, such neutron effectively scatters off the background nuclei if the mean-free 
time is shorter than the lifetime. On the contrary, if the mean-free time is longer than 
the lifetime, most of the energetic neutrons decay and become protons. Since the proton 
can be stopped in the thermal bath more easily than the neutrons, the hadrodissociation 
processes are suppressed in this case. 

Once the energy-loss processes via the electromagnetic interactions become ineffective, 
the effective lifetime of the neutron with energy En is given by 7„t„ (where 7„ = £^„/m„ 
is the Lorentz factor). Then, the condition for the mean- free time l/r^._^^^ being shorter 
than the effective lifetime is given by 




T»4.0xlO-keV(^) i::^:^) {^) , (5.7) 
and, for the cosmic time. 

Before closing this subsection, we check that the effects of the cosmic expansion are 
negligible; this can be seen if the cosmic expansion rate is smaller than the scattering rate 
of the nuclei. Indeed, in the radiation dominated universe with T^mg, the expansion 
rate is given by 



2 



^-(2-6sec)-x(^)\ (5.9) 
Thus, the condition T^'^^, ^ H is satisfied when 

T » 1.7 X 10- eV X /- m " f \ (5.10) 



Since we only consider the cases where X decays during the radiation-dominated epoch, 
this condition is automatically satisfied. 



5.3 Energy-loss of the hadrons via the electromagnetic interac- 
tions 

As we will discuss in the following sections, energy-loss of the energetic hadrons via the 
electromagnetic interactions is very important for the calculation of the abundances of 
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the light elements. If energetic hadrons completely lose their energy before scattering 
off the background nuclei, they do not cause the dissociations of the background nuclei 
(in particular, aBc)- In this case, only exothermic reactions are allowed for the hadronic 
particles, and the extra-ordinary hadronic inter-converting reactions between protons and 
neutrons become important. Such processes, which occurs at t < 100 sec, are discussed 
in Section El On the other hand, if the electromagnetic interactions do not stop the 
energetic hadrons (in particular, nucleons), they effectively scatter off the background 
nuclei and induce hadrodissociation. For t^lOO sec, this is the case. In particular, (i) 
ctBG is destructed and various energetic debris nuclei (n, p, D, T, and ^He) are produced, 
(ii) some of these energetic nuclei (in particular, T, ^He, and "^He) scatter off obg to 
produce heavier nuclei (^Li, ''Li, and ''Be), and (iii) the energetic nucleons emitted during 
the evolution of hadronic shower inter-convert the background proton and neutron through 
hadronic collisions. These processes will be discussed in detail in Section [7| and Section |H1 
In this subsection, we discuss the interaction of the energetic hadronic particles with 
the background (in particular, photon, electron, and positron) via the electromagnetic 
interaction. The energetic hadron Hi scatters off the background particles via the electro- 
magnetic interaction in particular with the following processes: the Coulomb scattering 
{Hi + —>■ Hi + e^), the Compton scattering {Hi + 7 — > Hi + '-/), the Bethe-Heitler 
scattering {Hi + 7 — ifj + e"*" + e~), and the photo-pion process {Hi + 7 — >• H'i+n). With 
these processes, energetic hadrons (gradually) lose their energy. The energy-loss rate of 
Hi can be expressed in the following form: 

\ dt J \ dt ) Coulomb WW CP WWbH \ dt /photo-pion' 

which is the sum of the various processes listed above. (The concrete expressions of the 
energy loss rates are given in Appendix El ) 

For our purpose, it is important to estimate how the hadrons lose their energy in 
the thermal bath before scattering off the background nuclei. If the energy-loss rate 
via the electromagnetic processes is large enough, hadrons are stopped before scattering 
off the background nuclei via the hadronic interactions. To estimate whether a hadron 
Hi is stopped or not through the electromagnetic interaction before scattering off the 
background proton or a, we calculate *^ 

l^-^A^Ei^^ E'n^T) ^ n^^ a%^^Jn. i^) ' dE,^, (5.12) 



and 



R""' {E^ ,K-,T)^T. (4f , -^T), (5.13) 



(depends also on the '*He fraction Y. We properly take into account the Y dependence in our 
numerical calculations. 
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where {dEnJit) is the energy- loss rate via the electromagnetic interactions which is given 
in Eq. ()5.11|) . and the sum in Eq. ()5.13p is over all the possible hadronic processes. Since 
p and a are the most abundant in the background among the nuclei, hadronic scattering 
processes are dominated by the scatterings with p and asc; in our study, we use the 
approximation 

^A,(T%^Ak = np(TH,+p^A^+... + no,aHi,+cc^A^, + -- (5.14) 

i 

R^^ {E^^^ , E'jj_ ; T) parameterizes the number of the hadronic scatterings before the en- 
ergy of the hadron Hi with its initial energy E'f^^ decreases to E'^. via the electromagnetic 
interactions. Thus, ii R"^{Ep^ ,E^^^]T)<1 (with eJ^^ being the threshold energy of 
some hadrodissociation process), the high energy hadron Hi is (mostly) stopped and does 
not cause the hadrodissociation process. On the contrary, if R^^ {E^f^\ E^^^^ ; T) is larger 
than unity, the high energy hadrons are not stopped through the electromagnetic interac- 
tion and cause hadrodissociation processes. In particular, if R^^{Ef\ E^^h T) ^ 1 with 
-^l/lo being the threshold energy for the destruction process of obg, high energy projectile 
hadron with its initial energy E'f^^ effectively destroys aBG- The number density of the 
obg becomes abundant after the cosmic time t ~ 200 sec (T ~ 0.1 MeV). Thus, if X 
decays after this epoch, abundances of the light elements (in particular, D and ^He) may 
be significantly changed by the direct destruction of aBQ. 

Using the quantity R^^ given in Eq. (|5.13|) . we estimate the energy of Hi just before 
scattering off the background proton or obg- For stable particles, such a energy for a 
given initial energy E^^^ , which is denoted as E^h~^\e^^^), is estimated by solving the 
following relation 

R"^{Ep^,E^^='^-T) = l. (5.15) 

For the neutron, we should take account of the fact that the neutron may decay before 
scattering off the background nuclei. Thus, for the neutron, we define E!^^^^\Ej^^^) as 
follows. For the given initial energy of the neutron we calculate the Lorentz factor 

7^™) = E^^ynin as well as the total scattering rate for the hadronic processes T^^^ = 
Upan+p^- + Uaan+a^-- (Here, the cross sections are estimated with the initial energy 
of the neutron.) Then, if > (in^^ Tn)'^ , we expect that the neutron scatters off the 
background nuclei before it decays, and we estimate E^^^='^\E^^''^) using Eq. ()5.15p . On 
the contrary, if < (7i"''Tn)~^, neutron is likely to decay before scattering off the 
background nuclei. In this case, the energetic neutron is equivalent to the proton for the 
calculation of the hadrodissociation processes. Thus, for this case we use the relation 

^(«=l)(^(in))^^(i^=l)(^(in))_#9 

^^Neutron loses its energy before it decays and hence the energy of the proton produced by the neutron is 
not exactly equal to En^\ However, the energy-loss rate of the neutron via the electromagnetic processes 
is much less efBcient than that of the proton and hence we can neglect the energy-loss of the neutron 
before it decays. 



27 




Figure 11: Cross sections of the pp scattering processes as functions of the kinetic energy 
of the beam (projectile) proton in the laboratory system. The solid line is the total cross 
section while the dashed (dotted) line is the elastic (inelastic) cross section. 




Figure 12: Cross sections of the np scattering processes as functions of the kinetic energy 
of the beam (projectile) proton in the laboratory system. The solid line is the total cross 
section while the dashed (dotted) line is the elastic (inelastic) cross section. 
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Figure 13: Cross sections of the pa scattering processes as functions of the kinetic energy 
of the beam (projectile) proton in the laboratory system. The solid line is the total cross 
section while the dashed (dotted) line is the elastic (inelastic) cross section. 
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Figure 14: Cross sections of the inelastic pa scattering processes. 
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For various nuclei, we calculate Ej^. as functions of the initial energy Ef^\ In 

order to precisely calculate E\j. , we need informations about the cross sections for the 
hadronic processes. For energetic proton and neutron, we use detailed experimental data 
of cross sections for pp, np, and pa collisions jTlj. In Figs. [TTJ [T21 [121 andHU we plot the 
data of total, elastic and inelastic cross sections for these collision processes, respectively. 

In Fig. E[ we plot the contours of the constant E^^^^^ on the T vs. -E^"-* plane. Here 
we use ?7 = 6.1 X 10^^° and Y = 0.25. From this figure, we see that the protons are 
completely stopped when the temperature is sufficiently high. This is because the number 
of the background electron is not Boltzmann-suppressed at high temperature, resulting in 
enhanced energy-loss rate. 

For a more quantitative discussion, it is convenient to define the typical energy-loss 
rate through the electromagnetic processes: 

If is larger than the hadronic scattering rate, energy loss is effective and the energy 
of the hadron is significantly reduced before scattering off the background nuclei. On the 
contrary, if F^m ~ ^A^^Af^y energy loss is not important. 

At the temperature T ~ (9(10 — 100) keV, Coulomb scattering is the most important 
for the energy loss. Using the energy- loss rates given in Appendix |B1 we can estimate T^u 
for the case where the Coulomb scattering is the dominant process: 

KMlcouiomb - (1-4 X 10-^Sec)-^ X /3^-iA (^) x;^/^^-..^ ^5 

where Xe = nie/T. (Here, we consider protons with (3p larger than the thermal velocity 
of the background electron, in which we are mostly interested at this temperature.) Since 
the number density of the electron becomes smaller as the temperature gets decreased, 
T^u in this case decreases as the temperature becomes lower. Comparing Eq. (|5.17p with 
Eq. (|5.2|) . we can see that the proton is completely stopped when C(10) keV. 

For ultra-relativistic protons, in fact, the Compton energy-loss (and the Bethe-Heitler 
process) becomes effective in particular at lower temperature. Typical behavior given in 
Fig. HSlis indeed understood by using the formula for, e.g., the Compton process. Using 
Eq. ()B.10|) . we can obtain F^m for the ultra-relativistic case: 

Km1c„.p» - (13 Bee)-' X (^) {^f . (5.18) 

Thus, for the ultra-relativistic protons, we can see that the energy loss is effective for the 
temperature T> C(O.l) keV x (^kin/100 GeV)-^ 

For the neutron, we plot the contours of constant E!^^^^ in Fig. ^| For the na 
scatterings, we use the cross section for pa processes assuming the isospin symmetry as 
well as small number of data for na scatterings. Empirically, the size of the isospin 
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Figure 15: Contours of the constant E^^^^\ the energy of the proton just before it scatters 
off the background nuclei {p or ^He). The horizontal axis is the temperature while the 
vertical one is the initial energy of the proton E^^^\ The numbers in the figure are the 
values of E^^^^\ and the solid line with "< -Eth.pa" shows the contour for E^^^^^ being 
equal to the threshold energy for the destruction process of '^He, -E"^^^ = 25 MeV. The 
horizontal thin-solid line denotes -E^"-* = E^^\ Here we use r] = 6.1 x 10~^° and Y = 0.25. 

breaking in this case is estimated to be within 20 % (10%) for E^""^ = 25 MeV (100 MeV). 
As we will discuss later, we will adopt 20 % errors to all the hadronic cross sections in 
our Monte Carlo analysis, which will also take account of this uncertainty related to the 
isospin breaking. 

As one can see, non-relativistic neutrons are effectively stopped when the temperature 
is so high that the number density of the background electron is large enough. For the 
complete stopping of the neutron, the temperature should be higher than ~ 100 keV, 
which is slightly higher than the stopping temperature of the proton; using Eq. ()B.19|) . 
one can calculate F^j^ and see that F^j^ for non-relativistic neutron becomes smaller than 
the hadronic scattering rate when T ^ 100 keV. Thus, at lower temperature, energy-loss 
of the neutron becomes inefficient. When the temperature becomes low enough, however, 
time scale for the hadronic scattering becomes longer than the lifetime of the neutron. 
This is the reason of the drastic change of E^^^^ at the temperature T ~ 0.3 keV for 
non-relativistic neutron. 

We also calculate the typical energy of D, T, ^He and ^He just before scattering off the 
background nuclei. We plot the contours of E^~^\ e!^~^\ Ei^~^^ and Ei^~^^ in Figs. [T7\ 
UHllTni andEHl respectively. (For readers' convenience, we also plot the cross sections for 
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Figure 16: Same as Fig. but for neutron. (The horizontal thin-sohd hne denotes 
= E^^^'^ accordingly.) 
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Figure 17: Same as Fig. ^1 but for D. The solid line with "< -Eth.o" shows the contour for 
E^) being equal to the threshold energy for the destruction process D + pbg — > 2p + n. 
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Figure 18: Same as Fig. UHl but for T. The solid line with "< -Eth,Tp" shows the contour 
for being equal to the threshold energy for the destruction process of the projectile 

T through T + pbg, while the thick solid line is the contour for Er^ being equal to the 
threshold energy for the process T + obg ^Li -l- n. 

the pD and p^He processes in Figs. |2I]and|22l respectively.) For these charged particles 
Ai, behavior of E\. is similar to the case of the proton. 

In this paper, we approximate that all the hadrons Hi have the energy Efj_ just 
before they scatter off the background proton or aBc- 

6 Inter- Conversion between Neutron and Proton 

So far, we have discussed the propagation of the hadronic particles in the expanding 
universe, paying particular attention to their energy-loss. Those particles cause various 
hadronic processes before and after being stopped by the electromagnetic scatterings and 
affect the abundances of the light elements. 

The first effect we discuss is the inter-converting p ^ n reaction caused by hadrons 
emitted from decaying X. Implications of such effect was first considered in Ref. [7]. 
In our study, we basically follow the strategy of Ref. [7] with several improvements, which 
are discussed in this section. 

Especially, for relatively early epoch T^O.l MeV (i.e., t^lOO sec), the emitted 
hadrons extraordinarily inter-convert the ambient protons and neutrons by the strong 

'^^"Of course, pions generated from the high-energy photons via the process 7 + 7bg ^ 7r+ + tt" are also 
included in the present work. 
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Figure 19: Same as Fig. UHl but for ^He. The solid line with "< i?th,3Hep" shows the 
contour for E^^^ being equal to the threshold energy for the destruction process of the 
projectile ^He through ^Hc + pbg? while the thick solid line is the contour for Ei^~^^ being 
equal to the threshold energy for the process ^He + aBG ^Li -l- P- 

interaction even after the conventional freeze-out time of the neutron in the SBBN. Thus, 
for massive particle X with relatively short lifetime {tx ~ (10~^ — 100) sec), the inter- 
conversion effect may induce a significant change of the light-element abundances. In par- 
ticular, since the proton is more abundant than the neutron, n/p ratio is enhanced if such 
inter-conversion is effective, resulting in overproduction of ^He. Thus, in this section, 
we discuss the inter-conversion processes. In particular, we present the thermally averaged 
cross sections for the relevant processes. In our analysis, for all the inter-conversion cross 
sections given in this section, we adopt 50 % uncertainties when we perform the Monte 
Carlo analysis. 

6.1 Cross sections of hadron-nucleon scattering 

First, let us summarize the cross sections for the relevant reactions. As we have discussed 
in the previous sections, once a high energy hadron is injected into the thermal bath at the 
early stage of the BBN (more precisely, when 0.1 MeV), energy- loss processes via the 
electromagnetic interactions are very effective (except for the neutral kaon). Thus, even if 
the non-thermally produced hadrons (in particular, 7r^,K^, p, and n) are quite energetic 

^^■^D/H and ^Li/H are also sensitive to n/p even at t > 100 sec. This effect is studied together with 
hadrodissociation effect in the next section. 
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Figure 20: Same as Fig. [T21 but for ^He. The solid line with "< -E'th,4Hep" shows the 
contour for -EI^"^'' being equal to the threshold energy for the destruction process of the 
projectile ^He through ^Hc + pbg, while the thick solid line denotes the contour for Ei^~^^ 
being equal to the threshold energy for the process ^He + obg ~^ ''Li + ■ ■ ■• 




Figure 21: Cross sections of the pD scattering processes. 
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Figure 22: Cross sections of the p^He scattering processes. 

when they are produced, they are quickly thermahzed. In this case, the typical kinetic 
energy of these hadrons becomes ~ T. Thus, for the inter-conversion processes, only the 
exothermic reactions are relevant since the kinetic energy of the thermalized hadrons are 
expected to be too small to induce endothermic reactions. In addition, the inter-conversion 
processes occur mostly with very small velocity. 

Concerning the exothermic reactions, it is well-known that the cross section a is ap- 
proximately inversely proportional to the velocity f3 when the velocity is very small; namely 

becomes almost constant. Thus, we use the threshold cross section. The cross sections 
given in the following are from Ref. 7j. 

The thermally-averaged cross sections for the processes tt'^ + n ^ p + ir^ and tt^ + p ^ 
n + 71^ are given by 

{aP)ll^ = 1.7 mb, (6.1) 
{a/3);Zn = 1-5C^(T) mb, (6.2) 

respectively. Here, C]j^{T) is the Coulomb correction factor which takes account of the 
modification of the wave-function of the charged particle due to the Coulomb field. For 
the opposite-sign charged particles. Coulomb correction factor is given by [ 75] 
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where aem is the fine structure constant and /i^^ is the reduced mass of the hadron Hi and 
the target nucleon. 
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Next, we consider the effects of the nucleons and anti-nucleons directly produced by 
the decay of X. In our study, we only consider the case where the numbers of such p and 
n produced by the hadronic decay of X are individually the same as those of p and n. 

When the temperature is high enough, nucleons are stopped in the thermal bath by 
the electromagnetic interactions. (Such temperature is given by T > 20 keV and 100 keV 
for (anti-) proton and (anti-) neutron, respectively.) In such a case, we treat nucleon- 
anti-nucleon pair NiNi like a meson, following Ref. 0. Then the NiNi meson induces the 
inter-conversion processes: N^Ni + Nj Nj + ■ ■ ■. If the nucleons are thermalized, we can 
use the threshold cross sections: 



When the temperature becomes lower, on the contrary, nucleons cannot be easily 
stopped. Then, they also induce the hadrodissociation processes and hence effects of the 
anti-nucleons cannot be studied just by taking account of the conversion effects. Impor- 
tantly, however, once the hadrodissociation processes become effective, a large number 
of protons and neutrons are produced while the anti-nucleons are not produced in the 
hadronic shower. In this case, anti-nucleons directly produced by the decay of X become 
irrelevant since their numbers are much smaller than those of non-thermally produced 
nucleons. Thus, when most of the energetic nucleons scatter off the background nuclei 
before being stopped, inter-conversion by the anti-nucleons becomes unimportant. When 
the anti-nucleons are not stopped, it is difficult to estimate their energies with which 
the inter-conversion cross sections should be evaluated. Thus, we do not include the inter- 
conversion process due to the anti-nucleon in such a case. Notice that we have numerically 
checked that the resultant constraints on the properties of X does not change even if we 
include the inter-conversion by the anti-nucleons with the threshold cross sections at such 
a later epoch. 

Finally, we comment on our treatment of the kaons. Since the kaons have relatively 
long lifetimes, they may also contribute to the inter-conversion processes. Unfortunately, 
however, cross sections of many of the conversion processes related to the kaons are not 
available. In addition, the neutral kaon is hardly stopped in the thermal bath and 
hence it is difficult to estimate its effects on the inter-conversion without making some 
assumptions. In fact, in Refs. [71 E], effects of the inter-conversion by the kaons are 
studied with some assumptions and approximations. We have followed the method given 
in those references and estimated the effects of the kaons. Then, we have found that, 
with the procedure given in Refs. [ZllIS], our resultant constraint (i.e., upper bound on 
the parameter mxYx) from the over production of ^He becomes slightly severer. In our 
analysis, however, we neglect the inter-conversion effects of the kaons since the inclusion of 
the kaons introduce some uncontrollable uncertainties. Effects of the kaons are expected 




(6.4) 
(6.5) 
(6.6) 
(6.7) 
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to make the constraints severer and hence, in order to derive conservative constraints, our 
treatment of the kaons is justified. 

6.2 Formulation 

In this subsection, we formulate the time evolution equations with the late-time "meson" 
injection. As we have mentioned, the hadron injection at the beginning of the BBN epoch 
enhances the inter-converting reactions between neutron and proton, and the freeze-out 
value of n/p may be increased. The inter-conversion term in the Boltzmann equations for 
the number density of the nucleon N[= p, n) is written as 

= -BhTxnx J2 i^N-.N' - Kn'-^n) , (6.8) 

IC N' 

where K^^n' denotes the averaged number of the transition N ^ N' per one hadronic 
decay of X. 

The averaged number of the transition N N' for one hadronic decay of X is expressed 

as 

K^^^, = Y.N''^P^l^j^„ (6.9) 

where Hi runs over the hadronic species which are relevant to the nucleon inter-converting 
reactions (i.e.. Hi = 7r+, 7r~, nn and pp). In addition, N^' is the averaged number of 
Hi produced by one hadronic decay of X. Note that we assume that two hadronic jets 
are produced in one decay of X in this section. A^^' is shown in Fig. El as a function 
of 2i?jet, where Ejet is the energy of one hadronic jet. Furthermore, P^\,j^i denotes 
the probability that a hadron species H^ induces the nucleon transition N ^ N' and is 
represented by 

Pn'-^n' = t;^— -p^, (6.10) 

dec abs 

where F^^ is the decay rate of Hi, and T^^^ = Fj^^^ + T^^^ + T^^^ + Ffi,„ is the total 
absorption rate. 

7 Hadrodissociation of Background a 
7.1 Basic equations 

As was discussed in Section El when R^^ is more than unity, the hadronic scattering 
processes between the emitted high-energy nucleons and the background proton or "^He 

'^^^Here, we only consider the effects of the "mesons" (tt"*", tt", nn and pp) directly produced by the 
decay of X. Notice that, for the period when the inter-conversion becomes important, the background 
temperature is so high that the hadrons we are interested in are stopped before scattering off the back- 
ground nuclei. Thus, the hadronic shower does not occur in this case, and we do not have to consider the 
secondary production of the mesons and baryons. 
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Figure 23: Schematic picture of hadron shower induced by a high energy neutron (proton) 
which scatters off the background proton or the background a^c- 

(called aBc) become effective. In particular, obg can be destroyed and energetic nuclei, like 
D, T, ^He are produced. We call this type of the hadronic destruction "hadrodissociation." 

In order to study the effects of the hadrodissociation, we follow the approach of Ref. [S] 
with several improvements. In Fig. 1221 we show the schematic picture of the hadronic 
shower induced by a high-energy neutron and proton. Hereafter, we discuss how we study 
the processes shown in Fig. 1221 

Before going to the main discussion, however, let us comment on the treatment of our 
high-energy anti-nucleons. In our analysis, we neglect hadrodissociation processes induced 
by high-energy anti-protons and anti-neutrons since we do not have sufficient experimental 
data for the scattering processes between an energetic anti-nucleon and a nucleus. Such 
anti-nucleons may change the abundances of the light elements by dissociating background 
nuclei (as well as by the inter-conversion effect which was discussed in Section . We 
expect, however, that the resultant constraints do not change much by this approximation 
since the anti-nucleons are not produced secondarily in the hadronic shower. Of course, 
such energetic anti-nucleon directly produced by the decay of X may produce energetic 
hadrons by scattering off the background nuclei, which may evolve into a hadronic shower. 
Number of such a process is at most the same as that of the hadronic shower induced by 
the energetic p and n and hence the numbers of the hadrodissociation processes may be 
underestimated at most by a factor of two. 
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Although the hadronic shower contains multiple scatterings of the hadrons in the ther- 
mal bath, the evolution of the shower can be followed by taking account of two types of 
elementally processes, which are both discussed in the previous sections: one is the elec- 
tromagnetic processes through which hadrons gradually lose their kinetic energy, and the 
other is the hadronic scatterings which change the number of the hadronic particles. 

Let us see what happens to a hadron injected into the thermal bath with the initial 
energy E^^^ in more detail. Such a high energy hadron may be a direct decay product 
of X or may be a daughter particle produced in the hadronic showers. As mentioned in 
the previous section, once a high energy hadron is injected into the thermal bath, it loses 
its energy down to ~ -E'^"^^ defined in Eq. fl5.15|l . In our study, we approximate that 

the hadrons injected into the thermal bath with the initial energy E'f^^ scatters off the 
background nuclei (i.e., pbg and a-Qo) with the energy -Sj^"^^ given above. 

Since there are various hadronic processes, many possible final states exist even if we 
fix the initial-state particles. We specify the individual processes by denoting {i,j;k); 
here, i and j specify the projectile and target nuclei, respectively, while k is for the final 
state. Then, we approximate the probability at which the projectile hadron scatters 
off a background nuclei via the process k) as 



nj{T)a 



ii,j;k) 



J2^rniT)a(i^rn;l) 
ni=PBG.aBG I 



(7.1) 



where Uj is the number density of the background nucleus j, and cr(ij;k) is the cross section 
for the process k) as a function of E^^\ 

For each processes, we need to calculate the energy distribution of the final-state nuclei. 
The relevant final-state nuclei are n, p, D, T, ^He, and ^He in this case. If we consider 
the scattering process (i, j; k) in which projectile hadron Hi hits the target nucleon and 
produces the final-state hadron i/^, the energy distribution of is given by 

Jh^ [^h.,EhJ-9h, --7^ , {'■^) 

where Eh- here is the energy of Hi at the time of the scattering while Eh^ is the energy of 
Hk just after being produced by the scattering process. In addition, g^^i'^^ is the number 
of Hk produced in the process (z, j; k). (For more details, see Appendix Using fH^'''\ 
we can also obtain the total energy distribution of H^ after the hadronic scattering of Hi, 
which is given by 

GH.^HdEi^\EH,;T) = j:PlIh)iE^^T)fiif\E^'\EH,). (7.3) 

j=p,a k 



'^^^Notice that the probabiUty P^^j.j.^ depends not only on -E^""* and T but also on the *He mass fraction 
Y . Such F-dependence of P^j.j^-j is taken into account in the numerical calculations. 

'^-'^''in the hadronic scattering processes discussed so far, Li and Be are not produced. Treatment of those 
nuclei will be discussed in the next section. 
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Notice that, if the initial-state particle has some energy distribution FHi{EHi'-,T) (where 
Ep^ is the energy of Hi just after being produced), then the distribution of the final-state 
particle Hj^ is obtained as 

F'hSEh,:T) =J:[ dE^FHAl^^T)GH.^H,{E^.EH,-T). (7.4) 
i •' 

With the relations given above, we can recursively follow the evolution of the hadronic 
shower. To make some image, let us consider the hadronic shower induced by a primary 
energetic nucleon generated by the decay of X. (To the primary nucleon, we as- 
sign the generation number "0.") We denote the initial energy of the primary nucleon 
as -E''™o)- As mentioned, emitted nucleon loses its energy via the electromagnetic pro- 

cesses typically down to E^\^^ , which satisfies i?^'(£^^™o), -E'^^o)^^ 2^) — 1- Then, uf^ 

scatters of the background nuclei (pbg or ckbg) via the hadronic interactions. We call 
this "primary scattering." In our analysis, the primary scattering occurs with the energy 
^'■^1^ , and we obtain the energy distribution of the "1st generation" hadrons if^^"* as 

i 

G^(o)^^(i){E^^(l), Ejj(i);T). (Here and hereafter, the superscript for the hadronic species 

are to identify their generation number in the evolution of the hadronic shower.) Thus, 
once the energy distribution of the primary hadrons F^^^ are known, we can calculate the 
distribution of the first-generation hadrons by using the relation 

F^h]^Eh,;T) ^Y.j dE^ffFS^{E^i:^;T)GH,-.H,{E^S:\EHr,T). (7.5) 

In our study, the energy distributions of the primary (i.e., "0-th generation") hadrons 
are calculated by using the JETSET 7.4 Monte Carlo event generator. Similarly, the 
distribution functions for the l-th generation nuclei are recursively calculated by using the 
following relation: 

fSI{Eh,;T) dEHA~'\EH,;T)GH,-.H,{EH,,EH,;T), (7.6) 

i 

where F^j. is the distribution function of in Z-th generation. After a large number of 
multiple scatterings, energy distribution function of any hadrons for energy region above 
the threshold energies of the hadrodissociation processes becomes negligibly small and the 
hadrodissociation processes stop. 

7.2 Approximations 

Evolution of the hadronic shower can be in principle understood by the recursive procedure 
discussed in the previous subsection. In the actual calculation, however, it is difficult to 
obtain the resultant distributions of the shower particles without any simplification and 
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approximation. One reason is that the number of hadrons contributing to the hadronic 
shower is enormous so it is difficult to track all the energy-loss processes of those hadrons. 
In addition, for some of the hadronic processes, experimental data for the cross sections 
are not available. Thus, in this subsection, we explain how we solve the basic equations 
given in the previous subsection with some simplifications and approximations. 

Our primary purpose is to obtain conservative constraints on the properties of X. 
By adopting reasonable experimental and theoretical errors, the SBBN scenario predicts 
abundances of the light elements consistent with the observations, as discussed in Section 
121 Thus, the non-standard processes usually make the theoretical constraints inconsistent 
with the observations. In particular, if hadronic scattering processes with the background 
obg occurs too much, D and/or ^He are overproduced. In addition, non-thermal processes, 
which will be discussed in the next section, may overproduce ^Li. 

Importantly, for some processes, we do not have enough experimental informations. For 
such cases, we adopt some approximation or assumption so that the numbers of D, ^He, 
and ^Li produced by the non-standard processes are minimized, resulting in conservative 
constraints. Thus, we should note that, for some case, the resultant abundances of D, ^He, 
and ^Li obtained from our calculations are their lower bounds. 

First of all, as we mentioned before, we simplify the treatment of the target (back- 
ground) nuclei by approximating that the energetic hadrons scatter off only the back- 
ground proton or obg- This is because most of the background nuclei are in the form 
of the proton or ^He. Of course, some of the energetic hadrons may scatter off other 
light elements in the background which may be destructed by those processes. If such 
processes are effective, however, production of those light elements from the dissociation 
of the background a^c is far more important since Obg is more abundant than other light 
elements (except Pbg)- Thus, our constraints is not affected by our approximation on the 
target particles of the hadronic processes. 

Second simplification is that, among various hadrons generated by hadronic scattering 
processes, only p and n are used as projectile nuclei in the next-step hadronic process 
(except for the non-thermal production processes of Li and Be, which will be discussed 
in the next section). This is because most of the nuclei produced in the shower processes 
are p or n, and also because we could not find sufficient experimental data for the cross 
sections for other nuclei. One might worry about the effects of the energetic ^He. However, 
energetic '^He is much rarer than p or n in the hadronic shower since the energy transfer 
to ^He in the elastic pctBC scattering process is fairly small. Notice also that the cross 
sections for the inelastic p + obg —>■]? + « + ■■■ reactions are relatively small. Thus, the 
energetic ^He has small effect on the evolution of the hadronic shower. 

The hadronic scattering processes considered in our study are listed in Tables El and |3J 
The experimental data of the cross sections are summarized in Refs. |7T| [73] . In addition, 
in our Monte Carlo analysis, we adopt 20 % errors to all the hadronic cross sections. 

In some case, we could not find sufficient experimental data and adopt some reasonable 
approximations or assumptions. In particular, the hadronic cross sections for the energy 
of the projectile higher than ~ 20 GeV cannot be found except for some pp and np 
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Process 



I = n 



I = p 



reaction type 



(i, pbg; 1) n+ pbg -^n + p p + pbg ^ p + p elastic 
(i, Pbg; 2) n + pec ^ n + p + n p + pbg p + p + it inelastic 
ih Pbgi 3) n + Pbg ^ n + n + jc p + pec ^ p + n + n inelastic 

Table 3: Hadronic processes with background proton pbg 



Process 



n 



p 



reaction type 



[z, a 
{i, a. 
{i, a: 
{i, a: 
(i, a 



n + aBG ^ n + a 
n + aBG ^ D + T 

n + aBG ^ 2r2 + '^He 
n + a—>-p + n + T 
n + aBG ^ n + 2'D 

n + aBG ^ p + 2n + D 
n + ttBG — ^ 2p + 3n 

n + ttBG n + a + 



V + «BG ^P + a 



P + aBG 
P + aBG 
p + a - 

P+OtBG 
P + "EG - 

P + aBG 
p + aBG - 



D + ^He 
p + n + ^He 
^2p + T 
^p + 2D 
'2p + n + D 
3p + 2n 
p + a + IT 



elastic 
inelastic 
inelastic 
inelastic 
inelastic 
inelastic 
inelastic 
inelastic 



Table 4: Hadronic processes with background aBG 



reactions. Fortunately, according to the existing data, however, the cross sections for the 
pp and np reactions are known to become approximately constant at high-energy region 
|71j . Thus, we assume that the inelastic cross section for the pa process is constant 
above E > 20 GeV. Our results are insensitive to this assumption since the evolution of 
the hadronic shower is mostly determined by the hadrons with energy less than ~ 0{1) 
GeV. In addition, sufficient experimental data are not available for the na reactions. For 
these processes, we use the SU(2) isospin symmetry and use the cross sections of pa 
reactions for the na reactions. Those cross section differ due to the Coulomb corrections. 
Using the familiar formula of the Coulomb correction factor [73], however, the Coulomb 
correction is estimated to be less than a few percent for the projectile energy larger than 
the threshold energy for the inelastic na scattering process (~ 25 MeV). Thus, we neglect 
such a Coulomb correction. 

In addition, the experimental data of the hadronic scattering processes for other pro- 
cesses are also insufficient. Thus, we make the following approximations for the daughter 
nuclei D, T, ^He, and ''He. 

(i) In considering the hadronic process, the energetic daughter nuclei scatter off only 
the background proton and aBG- 

(ii) The daughter nucleus is assumed to survive only if (a) its typical energy just 
before the hadronic scattering (i.e., ) is smaller than the threshold energy for 
the dissociation of by scattering off the background proton, and (b) typical energy 
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of the background photon in the center-of-mass frame is smaller than the threshold 
energy for the photodissociation process of Aj^. (In fact, the second condition is 
not crucial; the resultant constraints on X do not change much even if we do not 
include the condition (b).) That is, the surviving probability of the daughter particle 
Ak{= D, T, ^He or '^He) is simply given by 



[ 1 : otherwise 

where -E^^'^"* is the threshold energy for the hadrodissociation process, while Qa^. is 
the binding energy of A^. 

(iii) The daughter nucleus A^ is completely destroyed into energetic nucleons if -Pa^^a 
0. (For e.g., "^He + pbg — * 3]? + 2n.) 



k 



The approximation (i) is quite reasonable since almost all the baryons in the universe 
at the epoch we are interested are in the form of the proton or obg- The assumptions 
(ii) and (iii) are justified since our purpose is to obtain conservative constraints; indeed, 
the numbers of the non-thermally produced D, T, '^He and ^He are underestimated with 
these approximations since the energetic nucleons produced by the dissociation of the light 
elements rarely dissociate obg to produce other light elements. 

Since the target particle is pbg or q;bg, we can rewrite P{i].k) given in Eq. ()7.1|) as 



where cr^°p^ and crj^°*^ are total cross sections for the Nip-QG and Nia-QG scattering processes, 
respectively. In our numerical calculations, distributions of the final-state hadrons (in 
particular, those of the light-elements) are calculated with this formula. 

Since most of the final-state energetic hadrons are proton and neutron (as well as light 
mesons), we adopt an approximation such that only p and n are used as the initial-state 
energetic particles which cause hadrodissociations of q;bg- With such approximation, it is 
rather convenient to assign the generation numbers only to p and n by "integrating out" 
effects of other nuclei; using the distribution given in Eq. ()7.3j) . we define the distribution 
of the p and n after taking account of the effects of other light elements as 

Gn^n' {E'ff' , -Eat' ; T) = Gn^n' {E^n^ , E^r ; T) 

+ E /^^A,G^^A,(i?S'^\i?A,;T)P(t^^.._^,^...)(E^,;T) 



g^,^^--'^'--\Ef:'\E,,), (7.9) 

where and A^' correspond to p and n, and 5'^'°^ -^n + jg ^]-^g energy 

distribution of N' in the reaction A^ + ■ ■ ■ ^ N' + ■ ■ ■ with the energy of Ak being E)^^ 
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Although we have included the second term in Eq. ()7.9|) . which take account of the effects 
of the nucleons generated from the secondary destruction of the non-thermally produced 
light elements, it is much smaller than the first term. Indeed, for one hadronic decay of X, 
numbers of p and n produced in the shower processes are of (9(10 — 100) while the numbers 
of the destructed light elements are a few or less. Thus, our result is in fact insensitive to 
the approximation (iii) mentioned above. Furthermore, for other light elements Aj = D, 
T, ^He, and ^He, we define 

Gr,^A,{Ep\ Ea/,T) = Gr,^A,{Ep\EA^-T)PA^^A,{EA^;T). (7.10) 

Then, we define the distribution function of the light element Aj after the hadronic 
scattering of l-th generation nucleons, which we denote Fa ■ Notice that Fa obeys the 
recursion relation similar to Eq. ()7.6|) : 

F^}]{Ea^-T) = E JdENFt'\E^;T)Gr,-.A,{E^,EAr,T), (7.11) 

N=p,n 

where, in the above expression, Aj denotes all the possible light elements. 

In our numerical analysis, scattering and energy-loss processes of the energetic nuclei 
are studied by using Gn^ni and GAi^Aj defined above. (For details, see Appendix iDl) 

With the distribution function given above, we calculate the number of produced or 
destructed nuclei by the decay of X. In our calculation, target of the energetic hadrons 
is pbg or obg- In this case, numbers of the h = D, T and ^He are always increased by 
the decay of X (as far as we neglect the subsequent thermal processes). Then, with the 
distribution function, we calculate the total number of the nuclei Aj = D, T and ^He 
produced by the hadronic decay of one X: 

oo „ 

U,iT)=J2 dEA,Fj;]iEA,;T) 

= jdEA^jdEN E Sn{En-^T)Gn^aAEn,Aj-T). (7.12) 

N=n,p 

where 

oo 

SM{EN]T)=Y^F'ff{EM-T). (7.13) 

We also calculate the total number of destroyed ^He as 
UT)= E j dEN~SN{EN\T) 

N=n,p 

' ^ N f ~ 

E ^iN,a;k)iEN', T) + dEi^leG N^^UeiEN , Ei^le', T) {1 - P-i-Rc^iHciE^Ue'i T)} . 
.k=2 '' 

(7.14) 
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We found that the hadrodissociation of the '^He is dominated by the direct destruction 
in the hadronic process. Thus, although we have included the effects of the secondary 
destruction of the ^He which is made energetic via the first hadronic scattering processes 
(i.e., the second term in Eq. ()7.14|1 ). such effect is subdominant and does not significantly 
change the constraints. 

We should also consider the effects of the low-energy neutrons produced in the hadronic 
showers. Such neutrons may be energetic when they are produced, but they lose their en- 
ergy as they propagate in the thermal bath. (The energetic neutrons mainly scatter off the 
background p and ttBG-) Once the energy of the neutron becomes lower than the threshold 
energy of the destruction processes of the background obg (-E'na^ ~ 25 MeV), it no longer 
destruct the background obg- However, even after being thermalized, such extra-produced 
n may affect the abundances of the light elements. In particular, p may capture such low- 
energy n and D may be produced. Furthermore, ''Be may be dissociated by the thermal 
neutron via the process ^Be(n, ^He)'^He, which reduces the resultant abundance of ^Li. 
These processes are included in our BEN code. Importantly, at high enough temperature, 
the neutron with energy lower than E^^^ does not decay before being thermalized since 
its lifetime is much longer than the thermalization time. Thus, effects of the low-energy 
neutron is taken into account by injecting thermal neutron into the thermal bath. The 
number of the neutrons produced by the hadronic decay of X as 



In Figs. 1211123 and Uni we plot ^d, ^t, ^^he and as functions of the temperature T for 
rux = 100 GeV, 1 TeV and 10 TeV. As one can see, the ^-parameters almost vanishes at 
0.1 MeV. This is because, for such high temperature, energetic hadrons are stopped 
by the electromagnetic processes before scattering off the background nuclei. As the 
temperature becomes lower, the ^-parameters increases until T ~ O.SkeV. In this period, 
the hadrodissociation is dominated by the energetic neutron since the mean-free-path of 
the neutron is much longer than that of proton. Energy-loss of the neutron becomes less 
efficient as the temperature becomes lower, so the effects of the hadrodissociations become 
more effective as the temperature becomes lower. Once the cosmic temperature becomes 
lower than T ~ O.SkeV, however, the neutron decays before scattering off aeo- Since the 
stopping process of the proton is more efficient than that of the neutron, hadrodissociation 
is suppressed at the low enough temperature. Thus, we see sharp drop-off of the ^- 
parameters at T ~ O.SkeV. 

We also plot ^„ in Fig. 123 Note that, in Fig. 123 we subtract the number of neutrons 
which are contained in the initial spectrum of neutrons from ^„ in order to show the 
number of the secondary produced neutrons. The drastic decrease at T ~ O.SkeV is, again, 
due to the neutron decay. 




(7.15) 
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Figure 24: Number of produced or destructed hadrons per one hadronic decay of X with 
mx = 100 GeV. (solid line) is the number of the destructed ^He while (dotted line), 
(dashed hne), and ^SHe (long dashed line) are the number of D, T and ^He produced 
by X-decay, respectively. We take Yp = 0.25 rj = 6.1 x 10~^°, and 2Ejet = ^x- 




Figure 25: Same as Fig. |211 except for mx = 1 TeV. 
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Figure 26: Same as Fig. |21l except for mx = 10 TeV. 
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Figure 27: ^„ as a function of the temperature. Here, we take = 0.25 and r/ = 6.1 x 
the total energy of the two hadronic jets is 2Ejct = 100 GeV, 1 TeV and 10 TeV. 
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Figure 28: The cross sections of the ^Li production processes. The sohd (dashed) hne is 
for T + ^He ^ ''Li + n (^He + ^He ^ ^Li + p). 

8 Non-thermal Production of Lithium and BeryUium 

In this section, we discuss the non-thermal production processes of Li and Be. As we 
have discussed in the previous sections, energetic T, ^He and ^He can be produced by the 
hadronic or photodissociation processes with the background obg- Such energetic nuclei 
may scatter off the background obg again and produce other nuclei, in particular, ^Li, 
^Li, and ^Be. Although these collisions are not so frequent, they are important since the 
observations severely constrain the primordial abundances of ^Li and ^Li. 

First, we consider the non-thermal production of ^Li by the energetic T and '^He. In 
this case, energetic T and ^He are produced by 

f T ■ ■ • 

p{n)+aBG ^ S sjjg^ ... . (8-1) 

and these T and ^He scatter off the obg to produce ^Li: 

T + aBG^^Li + n, (8.2) 
=^He + aBG ^ ^Li + p. (8.3) 

Once the energetic nucleus = T or '^He is injected into the thermal bath, it loses its 
energy via the electromagnetic interactions by scattering off the background electron and 
photon while it also scatters off the background obg- With the energy-loss rate (dEAjdt) 
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j;iven in the previous section, number of ^Li produced by one Ai is 



(in) / 7 7-. \ -1 

A, 



j naaA,+a^a^^u+-{EAjpA,, (8.4) 

where E^^^ is the initial energy of Ai, (3Ai is the velocity of Ai, and Ea^~^^ is the typical 
energy of Ai just before its collision with background proton or obg (see Sectional). In 
addition, aAi+a-BG^^^i+--- cross section of the process (j8.2p or (j8.3p . (We plot the 

experimental data of them in Fig. |2H1 [ZHl EZl CS]-) Summing up the contributions of T 
and ^He, the number of ^Li produced by the above process per one hadronic decay of X 
is given by 



Ai=T,3He ° 

/"^^"^ ( dE^ \ ~^ 

y^(^li) J j ^aCrA,+aBG-«Li+...(^Aj/5A,^6Li^6Li, (8.5) 



where Pau^ai^i is the surviving rate of the non-thermally produced ^Li, f^^ is the 
cumulative energy-distribution function of energetic T and ^He during whole period of 
the evolution of the hadronic shower. For the energy distribution of T and ^He produced 
by the hadrodissociation processes, we use the experimental data [Hj (see Appendix 1^). 
Importantly, experimental result suggests that the energy distribution of the final-state T 
(and ^He) is almost independent of the energy of the initial state energetic neutron. Thus, 
we use the following formula for Z^. (with = T and ^He) 

JaA^aJ - -T^ • (8.6) 

Fitting formula for the differential cross section obtained from the experimental data, 
which is used in our analysis, is given in Eq. ()C.25|) . 

The non-thermally produced ^Li is energetic with their kinetic energies of (9(1 — 
10) MeV and might be destroyed by scattering off the background nuclei (in particu- 
lar, protons) before it is thermalized. To estimate the number of ^Li destroyed after the 
non-thermal production, we calculate the surviving probability Pen^eLi. (For details, see 
Appendix El) Then, for the cosmic temperature we are interested in (i.e., r< lOOkeV) we 
have found that the surviving rate is very close to 1. Thus, almost all the non-thermally 
produced ^Li survive until being thermalized. (Same is true for ^Li, ^Li, ^Be produced by 
the non-thermal processes with energetic ^He, which will be discussed below.) 

Using the cumulative energy-distribution function obtained by following the evolution 
of the hadronic shower, we calculate the ^6Li ^'^^ parameters for various background tem- 
peratures (as well as other cosmological and model parameters). The numerical result is 

"^^^In fact, non-thermally produced ^Li can be also destroyed by the process ^Li(pBGj ''He)'^He after being 
thermalized. Such an effect is taken in account in the standard code of the BBN calculation we used |5;-i| . 
and hence is not included in PsLi^sLi. 
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Figure 29: ^6lj as a function of the temperature for 2Ejet = 100 GeV, 1 TeV and 
10 TeV. Here, we take Yp = 0.25 and 77 = 6.1 x 10"^°. 



shown in Fig. 1221 As one can see, ^e^i suppressed when T ^ 100 keV. This cut off is 

from the fact that, at such high temperature, all the hadronic particles are stopped by the 

electromagnetic processes before scattering off aeG- (In addition, at high temperature, 

surviving rate of ^Li is also suppressed.) In addition, we see a sharp drop-off of 

at T ~ 0.3 keV, which is due to the decay of the neutron during the propagation in the 

universe. 

Next we consider the non-thermal production of ^Li, ^Li, and ^Be through the collision 
of energetic ^He with background q;bg- Such energetic '^He is produced by the elastic and 
inelastic scatterings between high-energy nucleons and background a. The number of the 
non-thermally produced nuclei per one decaying X is expressed as 

^-^4He/4Hel-i^4He 



where, here, = ^Li, ^Li, and ''Be. In Fig. IHO the experimental data of the differential 
cross sections for ^He + ^He ^Li + p, ^Be + n, and ^Li + X are plotted [ZHl EOl EH EH • 
As will be shown in Appendix El the elastic scattering of high-energy nucleons with aeo 
(A^ + aBG N + a) is not important for these production processes of ^Li, ^Li, and ^Be. 
That is because the transfered energy to ^He in the elastic scattering is much smaller than 
the case of the inelastic scattering (A^ + obg N + a + n's), although the cross section 
is fairly large. 
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Figure 30: The cross sections of the ^Li, ^Li and ^Be production process. The sohd 
(dashed) hne is for ^He + ^He —>■ ^Li + p (^He + ^He — > ''Be + n), while thick sohd hne is 
the total cross section for the process ^He + ^He — > ®Li + ■ ■ ■. 

Although the expressions ()8.5p and ()8.7p have the same structure, estimation of ^a^^'' 
is rather difficult. This is because we only have insufficient data for the transfered energies 
to "^He in the non-elastic scattering processes to produce energetic "^He. In particular, for 
the process p/n + ^He p/n + ^He + vr + • • •, energy distribution of ^He in the final state 
is quite uncertain at around the threshold energy of this scattering process. The number 
of the non-thermally produced Li and Be, however, depend on the energy distribution 
of ^He. At the present stage, we have to conclude that the reliable estimation of the 
number of non-thermally produced Li and Be from the process ^He + aBG Li/Be + ■ ■ ■ is 
difficult. Thus, we will not include the non-thermally produced Li and Be from this class 
of processes when we derive the constraint on X. 

We can, however, estimate the number of Li and Be produced from the process ^He + 
obg Li/Be -|- ■ ■ ■ by adopting the energy distribution of *^He generated from our shower 
algorithm (see Appendix [n]): 

oo „ 

Ar^4He(i?Ar, -E'4He). (8.8) 

N=p,n 1=0 

Using this relation, we estimate ^7u'^\ and ^73^°^ The results are shown in Fig. 

EH The results indicate that, by comparing with Fig. 1221 the effects of "^He-aeG collision 
are less significant than those of the T-ttBc and ^He-ttBc collisions. In addition, we have 
checked that, even if we adopt ^6u'^\ ^7u'^\ and ^7^'^^ obtained above, no significant 
change of the resultant constraints on Yx is seen. 



52 



Logio(t/sec) 
2 4 6 8 10 12 



10- 



10- 



•w^ 10-5 



10- 



10- 



=111 

1 


1 1 1 1 1 1 1 


1 11 1 1 11 1 ^ 






^2E„^=10TeV 1 


_ 


i/\ Y — 

' ■ v -^-^ - - 

A t ^. 




1 x\ 


2B.^^=lTeV E 


Il 


1 1 


77 = 6.1x10-1°: 
7^ = 0.25 
1 1 1 1 



10- 



10- 



10- 



- 10- 



1 0.1 10-2 10-3 IQ-i lQ-5 IQ- 

Temperature (MeV) 



10- 



Figure 31: ^7j^'^\ and ^73^'^'' as functions of the temperature for 2Ejet = 100 GeV, 

1 TeV and 10 TeV. Here, we take Yp = 0.25 and 17 = 6.1 x 10 
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So far, we have discussed the case where the parent nucleus inducing the non-thermal 
production of Li and Be is produced by the hadronic scattering processes. Energetic nuclei 
are, in fact, also produced by the photodissociation processes of obg- Using energetic 
photons produced in the electromagnetic shower processes, energetic T and ^He can be 
produced by the processes 

jT + p 
^ + "^^^\3He + n • 

Then, using T and ^He produced by the above process, non-thermal production of ^Li is 
possible with the process ()8.2|) and ()8.3|) . 

We take account of these effects including the following term in the Boltzmann equa- 
tion: 



dt 



poo 

^*He L(th,T)_,^p(th,T) C?^7^*He(7,p)T(^7)/7(^7) 
7+aBG^T/3He+- +^^6 

^(th.T) ^Et 1 mHeaT(4He,n)6Li(^T)/3T 

POO 

+^*He / (th,3Ho)^,^(th,3Hc) dE^(^iHe{^^p)^He{E^)fj{Ej) 
J^ithMc) «£/3He I I '^4HeO"3He(4He,p)6Lil-£^3HeJP3He, 
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(8.9) 



where ' and Eq ' {E^ ' ^ and Eq ^ ) are threshold energies for the processes 
7 + ^He ^ T + p and T + ^He ^ ^Li + n (7 + ^He ^ ^Re + p and ^He + ^He ^ ^Li + p, 
respectively. In addition, crT(4He,n)6Li(-E'T) (c"3He(4He,p)6Li(-^3He)) is the cross section for 
the process T + ^He ^Li + n {^B.e + ^He ^Li + p) with i?T (-E'3He)being the energy of 
the injected T (^He). For these cross sections, we use the formula given in Ref. ^Hj; -E't 
dependence of o"T{4He,n)6Li(-ST) is shown in Fig. EHl 



9 General Results 
9.1 Outline 

In this section, we present our numerical results. In particular, we compare the theoretical 
predictions on the abundances of the light elements with the observations and derive 
constraints on the properties of X. 

In our analysis, we first calculate the evolution of the number density of X using 



dn 



X 



dt 



-SHnx - TxUx- 



(9.1) 



At each temperature, photodissociation rates are calculated by numerically integrating 
the photon spectrum and the relevant cross sections. In addition, we also calculate the 
^-parameters defined in the previous sections. 

Then, we obtain the Boltzmann equations for the light elements. Evolution of the 
nucleons N{= p,n) is governed by 



dn 



N 



dt 



dn 



N 



dt 



+ 



J SBBN 



dn 



N 



dt 



+ BhnxTx^N + 



- photodis 



dn 



N 



dt 



(9.2) 



ic 



For other light elements with atomic number 2 or 3 (i.e., Ai = D, T, ^He), we obtain 



dnAj 
dt 



dnA, 
dt 



+ 



SBBN 



dnA, 
dt 



+ BtnxTx^Ai, 



(9.3) 



photodis 



while for ^He, 



dm 



He 



dt 



dt 



+ 



SBBN 



dm 



He 



dt 



(9.4) 



photodis 



'^^^Effects of the non-thermal production of ^Li by T and ^He produced by the photodissociations become 
important at relatively low temperature where -E'en becomes smaller than the threshold energies. Thus, 
in Eq. H8.9|l . the threshold energies of the ^Li productions are used for the lower bound of the integration. 
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Figure 32: Contours of constant D/H on the tx vs. i^vis^x plane for mx = 1 TeV. Here 
we take Bh = 1, -Evis = ^x, and X is assume to decay into two hadronic jets with 



'jet 



mx- 



(D/H) 



SBBN 



Here, we take rj 
2.78 X 10"^ 



6.1 X 10 . In the SBBN, the theoreical predication is 



For ^Li, we include the non-thermal secondary production process discussed in the previous 
section and hence we obtain 



dn. 



dt 



dri' 



6Li 



dt 



+ 



SBBN 



dn 



6Li 



dt 



+ 



- photodis 



dn 



6 Li 



dt 



- 7+aBG^T/3He+-' 



(9.5) 



Here, the terms with the subscript "SBBN" represent the SBBN contributions to the 
Boltzmann equations (including the effect of the cosmic expansion) . 

In order to solve these equations, as we mentioned, we have modified the Kawano 
Code (Version 4.1, with the nuclear cross sections being updated), including the new sub- 
routines which take account of photodissociation, inter-conversion, and hadrodissociation 
processes. The photodissociation and hadrodissociation processes included in our analysis 
are summarized in Tables |21 El and IH In addition, our treatments of the inter-conversion 
and non-thermal production of Li are discussed in Sections El and |H1 respectively. 



9.2 Predicted light-element abundances 

To see how the abundances of the light elements behave, we estimated the light-element 
abundances using the center values of the cross sections and model parameters. In Figs. 
IS21 ESI 01 ESI and Eni we plot contours of D/H, ^He/D, Fp, ^Li/H, and ^Li/H, in (r^. 
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Figure 33: Contours of constant ^He/D. Cosmological and model parameters are the same 
as Fig. ESI In the SBBN, the theoreical predication is (^He/D)gg3j^ = 0.335. 




Figure 34: Contours of constant Yp. Cosmological and model parameters are the same as 
Fig- EH In the SBBN, the theoreical predication is (^)gBBN ~ 0.249. 
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Figure 35: Contours of constant ^Li/H. Cosmological and model parameters are the same 



as Fig. 1221 In the SBBN, the theoreical predication is (^Li/H) 
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Figure 36: Contours of constant ^Li/H. Cosmological and model parameters are the same 
as Fig. 1221 In the SBBN, the theoreical predication of the abundance is (^Li/H)gggp^ = 
3.81 X 10^1°. 
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^xYx) plane. In the caculations, effects of photodissociation and hadrodissociation are 
both included. In addition, we take mx = 1 TeV and = 1, and consider the case where 
X decays into two hadronic jets with with the energy 2Ejet = mx- 

As one can see, for Tx ^ 10=^ - 10^ sec. He abundance decreases as the primordial 
abundance of X becomes larger. This is because, as Yx increases, hadrodissociation and 
photodissociation processes of ^He become more effective and hence larger number of ^He is 
destroyed. In addition, since the destruction processes of ^He are followed by the creation 
processes of D, ^He, and ^Li, abundances of these light elements first enhanced as Yx 
increases. If Yx is extremely large, however, all the light elements are destroyed; in this 
case, abundances of D, '^He, and ^Li are also decreased. Contrary to D, '^He and ^Li, ^Li 
is hardly produced. Thus, for larger value of Yx, more ''Li is destructed and hence the 
abundance of ''Li decreases as the primordial abundance of X increases. 

On the other hand, for shorter lifetime (rx ^ 10^ sec), inter-conversion between the 
proton and the neutron becomes more effective. In this case, Y and D/H increase as Yx 
increases. 

In Fig. inni we can see a distinctive trend of decrease of ''Li/H at tx ~ 10^ sec and 
EvisYx ^ 10"^'^ GeV. That is because free neutrons produced by hadronic showers are 
captured by ''Be through ^Be(n, ^He)^He, which reduces the resultant abundance of ''Li. 
This phenomenon was also pointed in |17j . 

9.3 Constraints 

Now we show the constraints on the primordial abundances of X, taking into account the 
theoretical and observational errors. For the execution of Monte Carlo simulation in BEN 
computation, we should understand the error of all of the reaction rates concerning both 
radiative and hadronic decay processes. For radiative decay processes, we have shown 
them in Table El 

As for the reaction rates related with the hadrodissociation processes, which aredis- 
cussed in Section and Section |Hl we should estimate their errors in advance of executing 
the Monte Carlo simulation. In this work, we assume that the cross sections and the 
other model parameters of transfered energies, which are used for computing ^'s, obey the 
Gaussian distribution with their la errors. Computing ^'s with such errors sufficiently 
many times (in our caes, one thousand times), we statistically evaluate the errors of ^'s. 
Here we adopt 20% error for all the hadronic cross sections {6a/a = 0.2), which is larger 
than typical errors of their experimental data. For the errors of the transfered energies 
to nuclei in final states after the collisions (the inverse slope parameter Kt in inelastic 
nucleon-a scattering, and slope parameters Bs\ in elastic nucleon-nucleon and nucleon-a 
scatterings (see. Appendix [nj , we adopt 20 % errors. In addition, as we mentioned in 
Appendix O the ^ parameters do not change much with the variations of the inelasticities 
{np and Ka). Thus, we neglect their uncertainties. Furthermore, we have checked that ^'s 
are insensitive to the variation of rj. Thus, we also neglect its uncertainty in evaluating 
the errors of ^ parameters. 
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Figure 37: Estimated theoretical errors of the ^-parameters for mx = 1 TeV, Yp = 0.25 
and 7] = 6.1 X 10~^°. 

In Fig. EZIwe plot the errors of ^'s as functions of the temperature. These are the case 
of mx = 1 TeV, Yp = 0.25 and rj = 6.1 x 10^^°. We have checked that the errors do not 
change much even if we change the values of mx and Yp. Based on this result, we use the 
following errors in our Monte Carlo analysis: S^n/in = 0.15, 5^d/^d = 0.2, 6C,t/^t = 0.2, 

5e3He/e3He = 0.2, S^^a = 0.2, and ^^fff "'VdS''''^ = 0.3. 

For the hadron-nucleon inter-conversion reaction rate which was discussed in Section 
we adopt 50% error for each cross section because there are not any adequate experimental 
data for the uncertainties of cross sections. Therefore, we take the larger errors to get a 
conservative bound here, for the details, see Ref. |T5] . 

To study the abundance of the light elements with X, we have included the effects of 
the photodissociation, inter-conversion, and the hadrodissociation processes into the BBN 
calculation. In order to estimate the theoretical uncertainties, we performed the Monte 
Carlo simulation. Here, we follow the basic procedure explained in ^21- In addition, 
for the BBN calculation, we take account of the observational error of 1] reported by the 
WMAP collaborations |22 (see Eq. (IZTT)|l ). 

In our statistical analysis, with a given set of the model parameters, we calculate the 
theoretical values of the light-element abundances and calculate xf, the likelihood variable 
for individual statistical variable Xj. For Xj = {nD/nn), and Y, we use both the upper and 
lower bounds from the observation, and hence 

= (^th)2 ('^obs)2 = i^^/^u) and F, (9.6) 
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Figure 38: Upper bounds on mxYx at 95% C.L. for = 1 and mx = 100 GeV. The 
horizontal axis is the hfetime of X. Here, the hues with "D/H (low)" and "D/H (high)" 
are for the constraints ()2.1|) and ()2.2|) . respectively. The straight dashed line is the upper 
bound by the deviation from the Planck distribution of the CMB. 



where xf^ and are the center values of Xi determined from the theoretical calculation 
and observations, while af^ and 0"°^** are their errors, respectively. In our analysis, (cx*^)^ 
is calculated by the Monte Carlo analysis. Notice that the depends on the model 
parameters through xf^ and af^. For Xi = r^^2 ijiau/nYi) and logiQ[(n7Li/nH)] we only use 
the upper bound. In this case case, we define xf as 



f^th ™obs\2 



^ ' ^ ' ' forXi=r3,2, (neLi/rin) and logio[(?^7Li/nH)]. (9.7) 

: otherwise 



Notice that, contrary to the case of SBBN, we do not use the lower bound on (riTLi/nn). 
This is because we do not include the non-thermal ^Li production processes through a-a 
collisions. All the observational constraints on primordial abundances of the light elements 
have been summarized in Section |21 

In Figs. EH El and HOI we plot the resuhs of the analysis at 95 % C.L. (i.e., xl = 3.84 
for Xi = (no/nn) and Y; xl = 2-71 for Xi = rg^a, {neu/nn) and logio[(n7Li/nH)]) on the 
Tx VS. E^isYx plane for mx = 100 GeV, 1 TeV, and 10 TeV, respectively. Here, the 
hadronic branching ratio is unity, and X decays into two hadronic jets with the energy 
2Ejet = iTT-x- As mentioned in Section |21 the constraint with use of the highest observed 
value of D/H (Eq. ()2.2p ) is shown together with that obtained by taking our standard 
value (Eq. (j2.1|) ). One can see that the constraint from D/H changes by a factor 2 — 3 by 
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Figure 39: Same as Fig. EHl but for nix = 1 TeV. 




Figure 40: Same as Fig. EHl but for nix = 10 TeV. 
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different adoption of the observed value. We also plot the upper bound not to disturb the 
Planck distribution of the cosmic microwave background. As one can see, constraints on 
the combination E^isYx is quite insensitive to the mass of X. This fact imphes that the 
constraints for the case of 2Ejet 7^ can be roughly estimated by rescaling the bounds 
given in the figures. 

In order to see the dependence on hadronic branching ratio, we also show the results 
for Bfi = 10~^ and 5/j = in Figs. I^and l^ respectively. (Here, we considered only 
our standard value of the observed D/H (Eq. ()2.1|l ).) As one can see, with larger value of 
the hadronic branching ratio, upper bounds on i^vis^x become severer. In addition, even 
with a relatively small value of the hadronic branching ratio (i.e., = 10"'^), the hadronic 
decay mode may significantly affect the light-element abundances. In addition, in order 
to separate out the effects of the hadronic decay, we also calculated the light element 
abundance only taking account of the inter-conversion and hadrodissociation processes. 
The resultant constraint is shown in Fig. 021 Notice that, in this figure, effects of the 
photodissociation is not included so the situation is unrealistic; this figure is shown just 
for demonstration. 

'^^^Our constraint for the case of Bh = is in a reasonable agreement with the results obtained by previous 
studies (in particular, by Cyburt et al. IF"). Cyburt et al. did not consider the constraint on nsjjc/^-D 
and, for riD/f^H, observational constraint milder than ours is used (1.3 X 10-5 < (no/nH)"''" < 5.3 x IQ-^). 
If we adopt {fiD/nii)"^^ used in 16 , the difference between the upper bound on Yx from our analysis and 
theirs is within the factor of 3 or so and is quite mild. It may be due to the difference of the photon 
spectrum used in the analysis. (For the comparison, notice that we have normalized the yield variable Yx 
by the entropy density s, and also that i?vis — ■^^x is assumed in 16 .) 
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Figure 42: Same as Fig. EHl but for mx = 1 TeV and = (no hadronic decay mode). 
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Figure 43: Same as Fig. EH but for mx = 1 TeV and no photo-dissociation. 
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With the hadronic decay modes, the most significant constraint on i^vis^x depends on 
the hfetime tx'- 

• For 10~^ sec^ tx ^ 10^ sec, the inter-conversion processes are efficient. In this case, 
significant amount of p may be converted to n and, consequently, Y may be enhanced. 
In this case, the constraint from the overproduction of ^He is the most significant. 

• For 10^ sec ^ ^ 10^ sec, energetic hadrons (in particular, neutron) is hardly 
stopped by the electromagnetic processes and hence the hadro dissociation processes 
become the most efficient. In particular, in this case, non-thermal productions of D 
and ®Li provide the most stringent constraint. 

• For 10^ sec ^ rx ^ 10^^ sec, energetic neutron is likely to decay before scattering off 
the background nuclei. In this case, effects of the hadronic decay modes become less 
significant compared to the case with shorter lifetime. Furthermore, in this case, 
effects of the photodissociation becomes comparable to or more significant than the 
hadrodissociation. Then, the strongest constraint is from the overproduction of ^He. 

Importantly, for the case with relatively short lifetime (i.e., rx ^ 10'' sec), the hadrodis- 
sociation and the inter-conversion processes are the most important. Consequently, the 
constraints strongly depend on the hadronic branching ratio. In this case, the upper 
bounds on E^igYx are approximately proportional to By^. On the contrary, for longer life- 
time (i.e., Tx ^ 10^ sec), the most significant constraint is from the overproduction of ^He 
by the photodissociation of aec- Then, the upper bound on E^isYx becomes insensitive 
to Sh- 

Figs. EH ~ 1121 are our main results and one can read off the constraints on the primordial 
abundance of long-lived exotic particles. Since our analysis does not assume any special 
properties of X, our results can be applied to various classes of the long-live particles. 
In the next section, we will discuss one important application of our results, i.e., the 
application to the gravitino problem. 

10 Application to Unstable Gravitino 

In the previous section, we have derived constraints on the primordial abundance of the 
late-decaying particle X. Now, we apply our results to one of the most important cases, 
i.e., the case with unstable gravitino. In supergravity theory, the gravitino, which is the 
superpartner of the graviton, exists. The gravitino acquires a mass 1713/2 from the effect of 
the supersymmetry breaking. In large class of the models, the gravitino mass is comparable 
to or larger than ~ (9(100) GeV. Importantly, the interaction of the gravitino is suppressed 
by inverse powers of the (reduced) Planck scale and hence its interaction is very weak. 
Thus, if the gravitino is unstable, its lifetime becomes very long. This fact means that, 
if the gravitino dominantly decays into visible-sector particle and its superpartner, the 
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decay of the gravitino in the early universe may significantly change the predictions of the 
SBBN unless the primordial abundance of the gravitino is small enough . 

With the inflation in the early stage of the universe, the primordial gravitino is once 
diluted but it is produced after the reheating starts. Thus, even in the inflationary models, 
we may still have the gravitino problem We can calculate the abundance of the 

gravitino 13/2 as a function of the reheating temperature Tr, which is defined in this 
paper as 



\ 1/4 

Tk^I-^^^'A] , (10.1) 



10 



with Finf being the decay rate of the inflaton and is the effective number of massless 
degrees of freedom. (Here, we use g^, = 228.75.) In our analysis, we reanalyzed the 
gravitino production processes in the early universe. In particular, we have used the 
thermally-averaged gravitino production cross section given in Ref. |^, which properly 
takes account of the effect of the thermal mass of the gauge bosons, and numerically solved 
the Boltzmann equation for the gravitino production. (Thus, the following fitting formula 
is different from the previous ones given in Refs. [HI IH21-) The details are discussed in 
Appendix |Fl Then, we find that the gravitino abundance after the inflation is well-fitted 
by the following formula: 



Y3/2 ~ 1.9 X 10-^^ 



X 



T, 



R 



1010 GeV 



0.045 In 



T, 



R 



10^0 GeV 



0.028 In 



T, 



R 



1010 GeV 



(10.2) 



Importantly, the primordial abundance of the gravitino 13/2 is approximately proportional 
to the reheating temperature Tr. Consequently, for the inflation models with high reheat- 
ing temperature, the gravitino abundance may become so large that the light-element 
abundances are too much affected by the decay of the gravitino. Thus, we obtain a upper 
bound on the reheating temperature after inflation. 

In order to derive the upper bound on Tr, we have to specify the decay mode of the 
gravitino to calculate its lifetime, Bh, and so on. In this paper, we consider two typical 
cases. The first case is that the gravitino can directly decay into the pair of colored 
particles. In particular, we consider the case where the gravitino can decay into gluon 
and gluino pair: ip^ ^ g + g (see Fig. HJ, producing one hadronic jet with Ejet = |'^3/2- 
Assuming that this is the dominant decay mode, the lifetime of the gravitino is given by 

T-3/2(^^ ^ ^ + ^) ^ 6 X 10^ sec X [j^^) ' , (10.3) 

where we have neglected other decay modes (in particular, decay into the lightest neu- 
tralino and its superpartner). In this case, the hadronic branching ratio is expected to be 
very large, so we use Bh = 1. In our analysis, we calculated the hadrodissociation rates 
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with the approximation that the numbers of hadrons (i.e., proton, neutron, pion, and so 
on) from the single gluon jet are the same as those from the single quark jet. In addition, 
since the hadronization processes of the decay products of the gluino are uncertain, we 
neglected the effects of the gluino. Such a treatment of the gluino might have made the 
upper bound on Tr less stringent by the factor of ~ 2 or so. For the photodissociation 
rates, since the gluino is a colored particle, most of its initial energy is expected to be 
converted to that of radiation. Thus, even though it is expected that some fraction of the 
initial energy of gluino is carried away by the lightest neutralino (which is assumed to be 
the LSP), we use E'vis = ""^3/2 in this case. 

Direct decay of the gravitino into the colored superparticles may be, however, kine- 
matically blocked. Thus, we also consider the case where the gravitino dominantly decays 
into the photon and the lightest neutralino (which we call "photino" in this section): 

— i> 7 + 7. In this case, the lifetime of the gravitino is obtained as 

r3/2(V^M ^ 7 + 7) ^ 4 X 10« sec (y^^^) ' • (10-4) 

Even if the gravitino dominantly decays into a photon and a photino, the hadronic branch- 
ing ratio is non- vanishing since the quark-anti-quark pair can be attached at the end of the 
virtual photon line. (See Fig. 0) In this case, Bh is expected to be of order C(aem/47r), 
and hence we adopt Bh = 10~^. In addition, because almost half of the initial energy 
is carried away by the photino in this case, we use the relation ii^vis = ^1^3/2- Here we 
assume that the decaying gravitino produces two hadronic jets with Ejet = |^3/2- 

In Figs.^^andE21 we plot the upper bound on the reheating temperature for the cases 
where the gravitino dominantly decays into gluon-gluino and photon-photino pairs at 95% 
C.L., respectively. For comparison, we also calculated the light-element abundances for 
the case neglecting the effects of the hadronic decay modes (i.e., for B^ = 0), which only 
takes into account the effects of the photodissociation. Excluded region for such a case is 
also shown in the figures by the shaded region. 

If the hadronic branching ratio is large (i.e, Bh ~ 1), the constraint on the reheating 
temperature is very stringent. For example, in the gravity- mediated supersymmetry break- 
ing scenario, the gravitino mass becomes comparable to the masses of the superpartners 
of the standard-model particles, and hence 7/13/2 ~ C(IOO) GeV is expected. In this case, 
the reheating temperature should be lower than 10^ — 10'' GeV. In the anomaly-mediated 
supersymmetry breaking scenario |5H^, however, the gravitino mass can be two to three 
orders of magnitude larger than the masses of the squarks and sleptons. In this case, the 
upper bound is slightly relaxed, and we obtain Tr ^ 10'' — 10^° GeV. 

For the case where the gravitino dominantly decays into the photon-photino pair, the 
upper bound becomes higher. However, the constraint is still much more stringent than 
the case with Bh = 0. For 1713/2 ~ (9(100) GeV, the upper bound is given by 10^ — 10® GeV. 
With larger gravitino mass, the constraint becomes less stringent. Since we are assuming 
Bh = 10~^ here, however, colored superparticles should be heavier than 7713/2- Such a mass 
spectrum looks quite unnatural from the naturalness point of view. 
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Figure 44: Upper bounds on the reheating temperature as a function of the gravitino mass 
for the case where the gravitino dominantly decays into gluon-gluino pair. Here, we take 
Bh = I, -Evis = ^3/2, and Ej^t = ^1^3/2- The shaded region is the excluded region for the 
case with B^^ = 0. 



Note that the above upper bounds on the reheating temperature can be alleviated if 
the primordial gravitinos are diluted by the decay of particles (such as the thermal inflaton 
or the moduli fields). Although again gravitino may be also produced directly by the decay 
or after the reheating via the scattering processes, there are a set of model parameters to 
make the bounds milder and resolve the gravitino problem. (For the recent study, see, for 
e.g., Ref. [84 ). 

The upper bound on the reheating temperature provides significant information about 
the evolution of the universe, in particular, for the successful scenario of baryogenesis. 
Since the primary baryon asymmetry is diluted by the inflation, baryogenesis should occur 
at the temperature lower than Tr. For some scenario of baryogenesis, this becomes a very 
stringent constraint. In particular, for the Fukugita-Yanagida scenario for baryogenesis 
jHS] where the lepton asymmetry generated by the decay of the right-handed neutrino is 
converted to the baryon asymmetry via the spharelon process, mass of the right-handed 
neutrino should be heavier than 10^ — 10^° GeV jHS]- Thus, for successful scenario of 
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Figure 45: Upper bounds on the reheating temperature as a function of the gravitino mass 
for the case where the gravitino dominantly decays into photon-photino pair. Here, we 
take = 10~^, ii^vis = |"^3/2, and E'jet = |'^3/2- The shaded region is the excluded region 
for the case with _Bh = 0. 



leptogenesis, the reheating temperature is required to be higher than 10^ — 10^'' GeV. 
Based on our analysis of the BBN, however, such high reheating temperature is allowed 
in a very limited case, if the gravitino is unstable. Thus, in order to realize the Fukugita- 
Yanagida scenario, we may have to make some assumption on the property of the gravitino, 
like a very heavy gravitino, a light stable gravitino, or an invisible decay of the gravitino 
(like decay into the axion-axino pair). 

11 Conclusions and Discussion 

In this paper, we have studied the BBN scenario with exotic late-decaying particle X. 
If there exist an exotic particle with the lifetime longer than ~ 0.1 sec, it decays after 
the BBN starts. Decay products of X may cause various non-standard processes, like the 
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photodissociation and hadrodissociation of the background nuclei, and the inter-conversion 
between the proton and neutron. As a result, abundances of the light elements are affected 
by the decay of X. Since the theoretical predictions of the SBBN scenario is in a reasonable 
agreements with the observation, decay of such long-lived particle during/after the BBN 
may spoil the success of the SBBN scenario. 

We have derived upper bound on the primordial abundance X, paying particular at- 
tention to the hadronic decay mode of X; as we have seen, as the hadronic branching ratio 
i?h becomes larger, upper bound on the primordial abundance of X becomes smaller. In 
particular, for the case where the lifetime of X is 10~^ sec ^ rx ^ 10^ sec, inter-conversion 
induced by the emitted mesons and nucleons provides the most stringent constraint. For 
longer lifetime (10^ sec ^r^^ 10'' sec), hadrodissociation of the background a becomes 
effective, resulting in non-thermal productions of various light elements (in particular, D 
and ^Li). If we consider the case with the long enough lifetime (10'' sec ^ ^ 10^^ sec), 
^He is overproduced by the photodissociation and hadrodissociation of ctBG; which gives 
the most stringent constraint on the primordial abundances of X. (See figures in Section 

El) 

Since our analysis has been done without specifying the detailed properties of X, our 
results can be applied to various exotic long-lived particles. Some of the examples are the 
gravitinos, cosmological moduli fields, and the NLSP for the case where the gravitino is the 
LSP. As an application of our analysis, we considered the unstable gravitino and applied 
our results to the gravitino problem. In particular, for the cases where the gravitino dom- 
inantly decays into gluon-gluino and photon-photino pairs, we derived the upper bound 
on the reheating temperature after inflation. For the case of the gravity-mediated super- 
symmetry breaking, gravitino mass is expected to be ~ (9(100) GeV. In this case, even if 
the hadronic branching ratio is ~ O(10~^), the reheating temperature is constrained to be 
smaller than 10^ — 10^ GeV. If the gravitino mass is much larger than ~ 0(100) GeV, the 
constraint on Tr may be relaxed. With such gravitino, however, the hadronic branching 
ratio would be close to 1 since, in such a case, all the superpartners of the standard-model 
particles are expected to be lighter than the gravitino from the naturalness point of view. 
(Such a mass spectrum may be realized in the anomaly-mediated supersymmetry break- 
ing scenario.) For m3/2 ~ (9(10 — 100) TeV with i?h ~ 1, the upper bound is given by 
Tr < 10^ - 10^° GeV. (See figures in Section CUl) 
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A Photon Spectrum 

In this appendix, we summarize the properties of the photon spectrum generated by the 
high energy photon from radiative decay of a heavy unstable particle. (Details of the 
calculation of the photon spectrum is discussed in 0.) 

Once the high energy photon is injected into the thermal bath consisting of the photon, 
(non-relativistic) electron, and nuclei, the high energy photon induces cascade processes. 
We have calculated the photon spectrum taking account of effects of the following pro- 
cesses: 

• Injection of the high energy photon from the radiative decay of X 

• Double photon pair creation (7 -|- 7bg — > e+ -|- e~) 

• Photon-photon scattering (7 -|- 7bg —^7 + 7) 

• Compton scattering off thermal electron (7 -|- e^Q — 7 -|- e~) 

• Inverse Compton scattering off background photon (e"*^ -|- 7bg — > e^*" -|- 7) 

• Pair creation in background proton (and aec) (7 + Pbg ^ + + p) 

Here, the subscript "BG" indicates that the corresponding particles are in the thermal 
bath. In addition, in the process 7 -|- pbg —^e^ + e~+p, the background proton plays the 
role of the source of the electric field. 

The photon spectrum is determined by following the distribution functions of the pho- 
ton and electron, and fe- Since the expansion rate of the universe is much smaller than 
the scattering rates of the electromagnetic processes, the relevant Boltzmann equations to 
be solved are written in the following forms 
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(A.l) 



IC 



dt 



(A.2) 



J DE 



where the terms with the subscripts "DP," "PP," "PC," "CS," "IC" and "DE" denote 
the contributions of double-photon pair creation, photon-photon scattering, pair creation 
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Figure 46: The photon spectrum for the background temperature T = 10 eV, 1 keV, 
and 100 keV (from above). Here, E^^q = 100 GeV and the flux of the injected photon is 
[df^/dt]BE = 5{E^ - E^,o) GeV-^ while [dfJdt]nE = 0. 



in nuclei, Compton scattering, inverse Compton scattering, and the radiative decay of X, 
respectively. 

Importantly, for the energy region E^ <C E^^q (with E^^q being the energy of the in- 
jected photon from the radiative decay of X), the function f-yi^E^) is determined by the 
total amount of injected energy per unit time; with this quantity being fixed, f-yi^E^) is 
insensitive to £^7,0- Thus, in the BBN with the radiatively decaying particle X, resultant 
abundances of the light elements primarily depend on the lifetime of X and the combina- 
tion E^^qYx, but are insensitive to E^^q if the combination E^^qYx is fixed. (For the case 
with the hadronic decay processes, they also depend on the hadronic branching ratio.) 

In Fig. Uni we plot the photon spectrum f-yi^E^) for several values of the background 
temperature T. As we mentioned, first of all, f^{E^) is insensitive to the energy of the 
injected photon as far as the total amount of the injected energy per unit time is fixed. In 
addition, f^{E^) drastically changes its value at E.y ~ m\/22T. This can be understood 
as follows. Photons with high enough energy easily scatter off the background photon and 
create e^e~ pairs. Since the number density of the background photon is much larger than 
those of the electron and nuclei, the pair creation process is extremely efficient for photons 
with E^ ^ m1/22T the number density of the photons with such energy is suppressed. As 
E^ becomes lower, however, the process 7 -|- 7 — > e"*" + e~ is kinematically blocked, and 
hence f-y{E^) is no longer suppressed. 
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B Energy Loss Rates 



In this Appendix, we summarize the energy-loss rates of high-energy hadrons in the electro- 
magnetic plasma. Here, we consider the energy- loss of the hadron with mass mn = AniN 
and charge Z. In addition, -Ekin, 7, and /? are used for the kinetic energy, Lorentz factor, 
and velocity of the hadron, respectively. 

In discussing the energy-loss processes via the electromagnetic interactions, one of the 
most important processes is the scattering with background electron (and positron). In 
calculating the scattering rate, it is necessary to determine the number density of the 
background electron (and positron). In our analysis, for the sum of the number density of 
the electron and positron in the background (called rig), we use the following approximated 
formula which well agrees with the exact formula: 



Up 



3 

2-7 



m, 



T > rrif. 



j^s 3/2 



27T 



^-m,/T . rn^^T> me/26 



(B.l) 



1 - ^F^r^n^ 



T < me/26 



B.l Energy loss of charged particles 

First, we consider the energy-loss rates of charged particles (i.e, p, D, T, ^He, ^He, and 
TT^). For those charged particles, the Coulomb, Compton, and Bethe-Heitler scattering 
are important. 

First of all, we summarize the energy-loss rate through the Coulomb scattering. (For 
details, see Ref. |7j.) For the temperature higher than the electron mass (T^me), the 
energy-loss rate of relativistic charged particle is given by 



dE, 



(ch) 



dt 



TT 



- Coulomb 
COS 6'min)] 



T>me, E],in>mH, 



(B.2) 



where A = ln[2/(l — cos^min)] — 1 with ^min being the minimal scattering angle in the 
center-of-mass frame. In discussing the energy-loss process in the thermal plasma, ^min is 
determined by the requirement that the energy transfer to the electron in the comoving 
frame be smaller than the plasma frequency |lUlj which is given by 



Airap 



nip 



Thus, in our analysis, we use the approximated formula 



A ~ In 



(B.3) 



:B.4) 
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where /? is the required velocity to go to the center-of-mass frame and 7 = (1 — 
Notice that, when the incident charged particle is non-relativistic in the center of mass 
frame, 7/? ~ 7/3. For nonrelativistic charged particle, we obtain 



(ch) 



dt 



- Coulomb 



Z a^^T A 

9 mH 



(B.5) 



For lower temperature (T 



< 



nip 



the formula of the energy-loss rate changes. For 



relativistic charged particles, we obtain 
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T<me, E],in>mH. 



(B.6) 



For non-relativistic particles, we need special care since the energy-loss rate strongly de- 
pends on the relative velocity between the charged hadron and the background electron. 
Denoting the distribution function of the background electron (as a function of its velocity 
Pe) as fej^^iPe), the energy-loss rate of the non-relativistic hadrons is given by 
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:b.7) 



Notice that, the integral in the first term of the right-hand side of the above equation is 
the number density of the background electron with velocity f3e < P- Since the hadron and 
the electron are both non-relativistic, contribution of the background electron with jSe > (3 
more suppressed than the contribution of the electrons with /3e < 13. Then, neglecting the 
second term of the right-hand side of Eq. ()B.7|1 . we obtain 
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Coulomb 



me/5 



where 
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'TT Jo 



(B.8) 
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Another important process is the Compton scattering. In particular, ultra-relativistic 
charged particles may lose significant amount of energy by scattering off the background 
photons. For the Compton scattering process, the energy loss rate is approximately given 

by 
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We use this formula in our analysis. 

Next, we consider the Bethe-Heitler process (e.g., for the proton projectile p + 7 — 
p + e"*" + e~) |in2j . This process is important for relativistic charged nuclei with energy 
E>AmHme/e^ ~ l.GAGeV x (T/0.1 MeV)-\ where A is the mass number and is the 
photon energy. The fitting formula of the energy loss rate is given by 



dE, 



(ch) 
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2 ry2 4 
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:b.ii) 



where tq is the classical electron radius (= acm/^e), and /^^^ is the distribution function 
of the background photon which is given by 



7bg ' 



.TT/ exp(e/T) — 1 
In addition, the function in the integrand is fitted by |lU4j 
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where the coefficients are given by 
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ai = 0.8048, 02 = 0.1459, 03 = 1.137 x 10"^ 04 = -3.879 x 10 



60 = -86.07, bi = 50.96, 62 = -14.45, 63 



3' 



and 



ci = 2.910, C2 = 78.35, C3 = 1.837 x 10^ 



(B.14) 
(B.15) 

(B.16) 



Although we include the photo-pion process (p + 7 — ^ p{n) + tt) into our analysis, 
it is important only for highly-relativistic nucleons with energy E^ 771^171.,^ /e^ ~ 4.7 x 
lO^GeV X (T/0.1MeV)~^. The fitting formula of the energy loss rate is approximately 
given by |105j 
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where eo = mJl + 



2mA 



is the (approximated) threshold energy of the photon for the 



photo-pion process in the rest frame of nucleon, and a' ~ 6.8 x 10 ^^cm^/eV is a constant. 
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B.2 Energy loss of neutral particles 



Now, we consider the energy-loss of the neutral particles, in particular, the neutron which 
is the only relevant neutral particle for which the energy-loss rate should be discussed. For 
the neutron, scattering with the background electron with magnetic-moment interaction 
is the most important process for the energy-loss. (Even for the neutron, we call such a 
process as "Coulomb scattering. 



For T^rrie, we obtain 
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where gn — —1.913 is the neutron magnetic moment jTI]. Furthermore, for T^me, the 
energy-loss rate is given by 
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where Ge{xe) = Xg + + Gxg -|- 6 with Xg = rrie/T. 

As well as the case of charged particles, the photo-pion process (n + 7 — ^ n{p) + vr) is 
also important for highly- relativistic neutrons, see Eq. ()B.17|) . We include this photo-pion 
energy loss rate in our analysis. Note also that, because the energy loss rate through 
n + 7-^r;, + 7is smaller than the Coulomb energy-loss rate [3110], we do not include this 
process in estimating the energy loss rate of the beam neutron. 



C Energy Transfer 

In this appendix, we discuss the energy transfers into the scattered/produced particles 
after the elastic and inelastic hadronic scattering processes. There are three types of 
energetic nuclei after the processes; scattered target nucleus {p or a), beam nucleon after 



#i8-y^jjgjj the temperature become so low (i.e., r7ie/26) that the number of the electron in the comov- 
ing volume becomes constant, expression for the non-relativistic case given in Eq. HB.19() becomes invalid. 
For such a low temperature, the energy-loss rate of the non-relativistic neutron is given by 



dE„ 



Coulomb inmlf \meT 



dt 

At such low temperature, however, the energy-loss rate of the non-relativistic neutron is negligible. Thus, 



Eq. HB.19|I is enough for our analysis. 
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the scattering, and daughter nuclei which are produced by the destruction of a. In order 
to study the evolution of the hadronic shower, it is necessary to consider how the energy 
of the initial-state particles are transferred to these final-state ones. 

Throughout this appendix, = {Ei,p^) denotes the four-dimensional momentum of 
Ni (or Ai) with its mass rrii in the rest frame of the target particle (i.e., comoving frame 
of the expanding universe), and Ki = Ei — irti is its kinetic energy. The dashed quantities 
are those for the final state particles. In addition, the mass difference between the proton 
and neutron is not important here so we use the approximated relation rrtp ~ m„ = mjq. 
Furthermore, mass of a nucleus Ai is also approximated by the relation m^. ~ AA^mN 
where is the atomic number of A^. 



C.l Elastic scattering 

We consider the elastic scattering of an energetic nucleon Ni {p or n) with a background 
nucleus Aj (background proton or a) through the process Ni{pi)+Aj{pj) — + Ni{p'j)+Aj{p'j), 
where we denoted the four-momentum in the parenthesises. Then, pj = {mup 0, 0, 0) and 
p'j = {Ej,p'j), and the kinetic energy of the scattered target {p or a) is given by 

K; = E'^~mA^. (C.l) 

Using the Mandelstam variable t, which is given by 

t=ip'^-p,y = -2mA^K'^, (C.2) 

we obtain 



Kinematically, the maximal possible value of Kj is given by 



For the elastic scattering processes, the distribution of t is well approximated by the 
following form with the slope parameter Bgi 



a 



1 d(T(°l) 



1 c/a("i) 



dt 



exp(-5,i|t|), (C.5) 



where a("') is the cross section. 

For the process pp — *• pp, an accurate fitting formula of -Bsi as a function of the center- 
of-mass energy ^/s is given in Ref. f91 . The error of the fit is within 10 %. In Fig. we 
plot Bsi for pp ^ pp as a function of ^/s. We use the fitting formula given in Ref. [HI] to 
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Figure 47: Bgi as a function of the center-of-mass energy ^/s. The sohd (dashed) hne 
denotes the slope parameter for the pp pp {pp pp) process. Here, the fitting formula 
given in |Q1^ is used. 



calculate the slope parameter. We also adopt the same formula for the elastic processes 
in which the projectile is neutron (i.e., npec ~^ ^p)- Notice that the Coulomb correction 
factor for the cross section is within 6% for = 100 MeV (15% for Kn = 10 MeV). 
Therefore, the error of the transfered energy for np — > np should be within 20%. With 
Eq. ()C.5|1 . we can estimate the averaged transfered energy as 

For example, for the process ppbg PP with Kp = 5 GeV, we can use this formula to 
obtain (K'^) ~ 60 MeV. 

For the elastic pa scattering, we use Bsi = 28 GeV~^ [02], which is derived from the 
experiment with the energetic proton with its kinetic energy 695 — 991 MeV. Although we 
neglect the energy-dependence of Bgi for the elastic scattering processes with q;bG) such 
approximation would not significantly affect our results. This is because the averaged 
value of the transfered energy through pa pa process is {K'^ ~ 5 MeV and the energy 
transfer in this process is less efficient than ppbg PP- We also adopt this value in case 
of the elastic scattering of high energy neutron, na — > na. 

'^^^We use the fitting formula even for 2 GeV ^ \/s < 4 GeV although primarily the formula was derived 
for \/s^4 GeV |^. Note that the experimental data of the cross sections suggest that the energy 
dependence of Bsi is very mild. 
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Then, from Eq. ()C.5|) . we derive the energy distribution of the final-state Aj as 



where = 2mAjK'^^™^^^ is the maximal possible value of |t|. In our study, we use 

the above formula to evaluate the energy distribution of the scattered Ai. Notice that 
the energy distribution of the Ni in the final state can be derived by using the relation 

C.2 Inelastic scattering 

In this subsection, we discuss the inelastic scattering process of an energetic nucleon Ni 
with a background nucleus Aj. As we have discussed in the previous subsection, the 
differences of cross sections for the pp and np processes are small. Therefore we apply the 
formulae obtained for the pp and pa collisions to np and na collisions, respectively. 

Since we sometimes use the Lorentz transformation from the center-of-mass frame to 
the rest frame of the target particle, it is convenient to define the Lorentz factor as a 
function of the kinetic energy of the beam nucleon. Explicit expression of the Lorentz 
factor is given by 

(CM) , _ {mM^+mA^) + K, 

where the energy in the center-of-mass frame is given by 

E^^^){K,) = ^{mN,+mA^f + 2mA^K,. (C.IO) 

C.2.1 Inelastic scattering with background proton 

First, we consider the inelastic scattering between high-energy proton (neutron) with back- 
ground proton. The most important inelastic processes are the ones with one pion pro- 
duction: p + pbg P + p{n) + vr {oi n + pbg ^ n + p{n) + tt). Indeed, the cross sections 
of these processes are much larger than those of other inelastic processes. Hereafter, we 
consider the process p{pi) +Pbg(Pj) ^ p{p'i) + Nj{p'j) + 7r(p'^). 

In order to study the energy distribution of p in the final state (i.e., distribution of 
K-), it is useful for us to introduce Feynman's x-parameter. Let us decompose the three- 
momentum p. as 
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where y and p[^^ are components of p[ parallel to and perpendicular to the direction 
of the initial-state energetic proton, respectively. Then, the Feynman's x-parameter is 
defined as 

X = -^r^, (C.12) 

/(max) ' ^ ' 

where y = |p-y| and p'i^"^^^ is the maximal possible value of p-y. From experimental 
data on the high-energy inclusive pp reactions, it is known that the distribution of x is 
approximately independent of i/i (called "Feynman scaling" jHSHMl)- Indeed, various ex- 
periments reported that the kinetic energy of the final-state proton K'^ in the inclusive pion 
production process has flat distribution (see, for e.g., |HSj). Thus, with such experimental 
result, we estimate the averaged value of 

1 ^^(incl) 

. X — ax 

(x) = ^ — -4^^ 0.5. (C.13) 

; dx 

dx 

Using the above result, we parameterize the energy of the beam particle (i.e., p) in the 
final state using the "inelasticity" k^: 

Kl = {l-Kp)K,. (C.14) 

Notice that, for high-energy pion production processes, z ~ 1 — Kp. ^^'^ In our analysis, 
we approximate that Kp is constant and use Kp = 0.5 for the numerical calculations. 
That is, we approximate that the fraction of the energy-loss of the initial-state energetic 
proton is always Hp. We also use the above arguments in case of the energetic neutron 
injection. Then, we obtain the energy distribution of the final-state p as 

/Lbc;2)(^^' K'i) = /Lbg;3)(^^' ^0 = HK-a- ^p)K.) . (C.15) 

In order to study the energy distribution of other final-state particles, we can use the 
experimental data of the momentum distribution of the daughter particles from the pp 
collision. In particular, if we consider some daughter particle produced by the pp 
collision, energy distribution is well fitted by the following formula with the inverse slope 
parameter Kt jHSl IHZI : 

1 dcT^™''^) 

oc expi-MT/Kr), (C.16) 



MrdMrdydcf): 



^^'^ Although, strictly speaking, the relation x ~ 1 ~ Kp holds in the CM frame, it can be approximately 
satisfied even in laboratory system for high-energy collision. 

#^^We have checked that, for the 20 % variation of the K-parameters, ^-parameters change at most 15 
%. Such uncertainty is relatively small compared to the total uncertainty of the ^-parameters due to the 
errors of the hadronic cross sections (see section El. Thus, in our analysis, errors from the K-parameters 
are neglected. 
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where Mt is the transverse mass which is give by 

M^ = ml^ + \p'H^J\ (C.17) 

y is the rapidity and 0a is the azimuthal angle. Remarkably, for any kind of the daugh- 
ter particles such as protons, pions, and the other mesons, Eq. ()C.16|) provides a good 
approximation. In addition, the inverse slope parameter is not sensitive to ^/s. 

It is notable that Kt provides an typical kinetic energy of the daughter particles in 
the CM frame since (Mt) ~ Kt- This is interpreted as the consequence of the fact that 
the fragmentation scheme into hadrons is controlled by QCD and that represents the 
typical energy scale of the fragmentation. Using the results of the resent experiments of 
high-energy pp and p^Be collisions, we adopt the following value for Kt 



/sTt = 140 ± 15 MeV. (C.18) 

From the distribution of Mt, we estimate the distribution of the kinetic energy of the 
scattered target Nj. Unfortunately, direct information of is not available. Thus we 
assume that, at the high-energy scattering processes we are interested in, the daughter 
particles are distributed isotropically in CM system. With this assumption, we estimate 
the averaged value of the kinetic energy of Nj as 



{K'^),, = j^'^''\K,)Kt + [i^'^'^KK,) - Ij m^. (C.19) 

For simplicity, we approximate that the scattered target always have the kinetic energy 
{Kj)pp and hence the distribution of the scattered particle is written as 

/(™;2)(^^' = /(;PBc;3)(^^, = 5{K'^ - {K'^U). (C.20) 

In a similar fashion, we can also estimate the averaged kinetic energy of the produced 
pion as 



{K)pp = 7^^^n^.)^T + [7^^™n^.) - Ij (C.21) 

However, for pions, the time scale for the hadronic scattering processes is always longer 
than that of the electromagnetic stopping processes or the lifetime. Thus, we do not have 
to consider the energetic pions. 



C.2.2 Inelastic scattering with q;bg 

Next we consider the energy transfer in the inelastic processes with background obg- In 
this process, the target particle obg may or may not destroyed so we have to consider 
various final-states. As in the case of the inelastic scattering with the background proton, 
we denote the process as p{pi) + aBciPj) ~^ p{p'i) + ^jip'j) + '^(Pn) + " " "• 
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For the energy distribution of the injected particle in the final state, we use the similar 
treatment as the inelastic pp collision. That is, we define 

Kl={l~ K^)K,, (C.22) 

where Ha is the inelasticity related to ^He. As in the pp collision case, we approximate 
that constant and use the distribution function 

/rp,«BG;2)(^- ^0 = /rp,«BG;3)(^- ^0 = HK:-il- . (C.23) 

Since we do not have sufficient data to estimate Kq,, we assume that Kq, = Kp in our 
analysis. 

Now we consider the daughter nuclei produced by the inelastic scattering processes with 
a-QQ. In our study, energy distribution of the final-state nuclei is mostly approximated as 
follows. If kinematically allowed, we assume that all the produced debris nuclei (except 
the injected proton) have kinetic energy 



(ifOpa = 1^'-''KK,)Kt + [l^'-'^KK,) - Ij (C.24) 

where we used the transverse-mass distribution given in Eq. (HnH). (The total energy 
is assumed to be conserved by the pion emissions.) If Ki is not large enough, the total 
kinetic energy of the final-state particles becomes larger than Ki with the above approxi- 
mation. In such cases, we assume equipartition of the total momentum in the CM system. 
Therefore, our approximation is surely conservative than the case that we always assume 
the equipartition at every moment. 

In the study of the shower evolution and the surviving probability of the light elements 
produced by the hadrodissociation of q;bg! we use the energy distribution obtained by 
the above procedure. For some processes, however, experimental data provides better 
information about the energy distributions. In particular for the processes n + Obg ^ 
T H — ■ and n + obg '^He + ■ ■ ■, the daughter T and ^He play the role of the "spectator" 
particle and their typical energies are known to be very small (~ 0{1) MeV) [Zl]. For 
these processes, the energies of T and '^He are overestimated if we naively use our shower 
algorithm. For the process n + a — T + ■ ■ -, the experimental data is available and 
we found that the energy distribution is the final-state T is well approximated by 

1 rfa.^._T+... ^, Q Qgg3^ ^ , (C.25) 



a„+a^T+... dE^ V 5.789 MeV 



'^^^In the so-called "Cascade Model," where a nucleus such as ^He is approximately treated as a 
composition of many "independent" nucleons with uniform density, an approximate relation 
1 — (1 — KpY-^~^ '^^ holds (sec Ref. |1U(J| and references therein). If we adopt Kp ~ 0.5, we obtain 
Ka ~ 0.59 — 0.63 and the error in our treatment is about 20%. 

'^^^Since the energy of T and '^He may be overestimated in the study of the shower evolution, the energy 
of the final-state p and n may be underestimated in some case. This may be the case when the Ki is so 
small that the equipartition of the total momentum in the CM system is adopted. Then, such an error is 
expected to be small since, in such a case, energy of the nuclei are inversely proportional to their masses. 
In addition, our treatment is also justified from the point of view of obtaining conservative constraints. 
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independently of the energy of the beam particle. The energetic non-thermally produced 
T and ^He play very significant roles in the study of the non-thermally produced ^Li. In 
calculating the number of the non-thermally produced ^Li, we adopt Eq. (jC.25jl as energy 
distributions of T and ^He produced by the hadrodissociation processes. 

D Method for the Numerical Calculation 

In this appendix, we explain how we study the evolution of the hadronic shower in our 
numerical calculation. In our approximation, once a high-energy hadron Hi is injected into 
the thermal bath with the initial energy E'f^\ its energy is first decreased down to e')/^~^\ 
then scatters off the background nuclei (in our case, proton or a-Qo). Then, the energy 
distribution of the final-state particles after the hadronic scattering process is given by 
the distribution function Gh^^h defined in Eq. ()7.3|) (or Gat^^tv given in Eq. ()7.9|) with 
our approximation). The evolution of the hadronic shower can be in principle solved by 
using Eq. ()7.fij) once the initial spectra of the primordial hadrons directly emitted from X is 
given. Since our purpose is to count the number of the non-thermally produced/ dissociated 
light elements, however, there is rather a simple method which we use in our numerical 
calculation. In our numerical calculation, we first prepare the energy bins for kinetic 
energy of each relevant hadrons (in our method, we prepare the energy bins for proton, 
neutron and ^He). Hereafter, we consider the case where the J-th energy bin has the 
center value Ej and the width IS.E1. (Thus, Ej ± ^(Ai?/ + AEi±i) = Ej±i.) Here, / runs 
from to A^bin where 0-th bin is for stopped hadrons while iVbin-th bin is for hadrons with 
maximal possible energy. Then, we calculate the accumulated number of hadrons fallen 
into each energy bin during the evolution of the hadronic shower. 

As we mentioned in Section [7| in our approximation, evolution of the hadronic shower 
is studied by following the energy-loss processes of the proton and the neutron. Thus, the 
first thing to do is to calculate the number of the nucleon Nj scattered into J-th bin when 
one nucleon Ni is injected into J-th bin. By using G given in Eq. ()7.9|) . we can calculate 
such transfer matrix, which we call T^fjj^, as 

T^'n, = Gn^-.n,{Ei, Ej- T)AEj, (D.l) 

where E^]^~^^ satisfies the relation R^^{Ei, E^j^~^^ \T) = 1. Since, in our case, the energy 
of the nucleon always decreases after the hadronic scattering processes, T^fi^. vanishes if 

/ < J. #24 

Next, we consider the case where one Ni is injected into I-th bin, and estimate the 
accumulated number of Nj fallen into J-th bin after (infinite number of) multiple hadronic 
scattering processes. We call such quantity as UNi,Ny Using the fact that T^f^j^^ = for 

#^^Strictly speaking, j^, does not vanish since the energy bin has some finite width. In our analysis, 
such correction is taken into account. 
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I < J, the following relation holds 

i-i 

K=J+l Affe 

Thus, using the above relation, we can recursively determine Un^Nj from / = J = by 
increasing I and J. 

For the hadronic decay of X, then the accumulated number of Nj fallen into J-th bin 
is calculated as 

oo 

S^^iEj)AEj ^J2Fi?,m^Ej = Y.Y.Fn}{Ei)^EiU'/^^^^, (D.3) 

1=0 I Ni 

where F^'^ is the distribution function of Ni emitted from the hadronic decay of X. Sub- 
stituting this relation into Eqs. ()7.12|) and ()7.14|) . we calculate the ^-parameters for the 
light elements. 



E Dissociation of Non-thermally Produced Li 

Here, we discuss the scattering and dissociation of the non-thermally produced Li (and 
Be). Non-thermally produced ^Li and ^Li with energy of the order of (1 — 10) MeV can be 
destroyed by scattering off background particles before being stopped. Thus, in order to 
estimate the present abundance of ^Li and ''Li, it is necessary to understand the surviving 
probability of these nuclei. 

As other charged nuclei, the energetic Li loses its kinetic energy during the propagation 
in the thermal bath by scattering off the background e and 7. Since the scattering rate 
(and the energy-loss rate) of ^Li is much larger than the expansion rate of the universe, 
the evolution of the number density of the non-thermally produced ^Li is governed by 

1 (inLj(N.T.) 



(E.l) 



?T.Lj{N.T.) dt 

where au+p^--- is the cross section for the Li dissociation process. Thus, the surviving rate 
of Li with its initial energy E^^"* is calculated as 



Pi 



Li^Li 



{E 



(m)^ 
Li y 



exp 



exp 



dtUpCTu 



.(th) 




(3u 



(E.2) 



is the threshold energy for the 



where (dEu/dt) is the energy-loss rate of Li, and -^Li+p- 
destruction process of Li. 

#^^In fact, dissociation process ^Li(p, '^He)^He may occur even after ^Li is thcrmalized. Such an effect 
is taken into account in the BBN code we use and hence the surviving rate Pen^eLj parameterizes the 
dissociation of ®Li before thermahzation. 
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Figure 48: Cross section for the process ^Li(p, ^He)^He as a function of the initial energy 
of ^Li. For the original data, see Refs. [H7| I77j. 



In the calculation of the surviving probability of ^Li, the dominant dissociation process 
is ^Li(p, ^He)^He. Cross section for this process is given in Refs. [SZl HZ]. For the 
readers' convenience, we show the plot of the cross section for this process in Fig. 

We numerically evaluate the surviving rate for ^Li for various values of its initial 
energy. In Fig. 03 we plot the contours of constant surviving probability Pei^i^eu- As one 
can see, for T^50 keV when the thermal process ^Li(p, ^He)^He becomes inefficient, the 
surviving rate Pau^ei^^ becomes almost 1 for the typical value of the initial energy of ^Li 
(i.e., eI^^ ~ (5 - 10) MeV). 

For non-thermally produced ^Li, the cross sections for the processes ^Li(p, '^He)'*He 
and ^Li(p, ^He)X are available. We use the cross sections for these processes given in 
Refs. jHHl EZl EZI and Ref . jSfl [77] , respectively. Then, we calculate PrLi^TL; and found 
that PrLi^TLi is also close to 1 for the cases we are interested in. 

Finally, we comment on the surviving rate of non-thermally produced ''Be. Unfortu- 
nately, the cross sections for the dissociation processes of non-thermally produced ^Be is 
not available. Thus, we use the dissociation cross sections of ^Li also for the dissociation 
processes of ^Be. 

Although the surviving rate Peu^f^u is almost 1 for the case where the non-thermal 
production of ^Li may become significant, we include Pei^i^eu in the calculation of the 

'^^^There is another possible dissociation process of ^Li(p, 7)^Be. Cross section for this process is, however, 
much smaller than that for ^Li(p, '^He)'*He |77], and hence this process is irrelevant. 
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Figure 49: Contours of constant Peu^^^u- The horizontal axis is the background temper- 
ature while the vertical axis is the initial kinetic energy of ^Li. 



number of non-thermally produced ^Li. Same is true for ''Li and ''Be when non-thermal 
production of these light elements are discussed. 

F Primordial Abundance of Gravitino 

In this appendix, we evaluate the primordial abundance of the gravitino in the inflationary 
universe. In particular, we calculate the gravitino abundance as a function of the well- 
defined reheating temperature Tr. Compared to the old result given in Ref. 0, there 
are two major improvements: (i) we have included the effects of the gravitino production 
during the period when the infiaton field is still oscillating (i.e, in the so-called dilute 
plasma), and (ii) for the gravitino production cross section in the thermal bath, formula 
taking account of the thermal mass of the gauge bosons is used in order to avoid the 
infrared singularity. 

Gravitino is the gauge field for the local supersymmetry, and hence it couples to the 
supercurrent as 
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where xr and are fermion and boson in chiral multiplets, A is the gaugino, and F^p is 
the field strength of the gauge field. (Here, denotes the covariant derivative and xb. 
satisfies (1 - '^^)xr = 0.) 

Importantly, the interaction of the gravitino is suppressed by the inverse of M* and 
hence the gravitino interacts very weakly compared to the standard-model particles. In 
particular, if the temperature is lower than ~ M*, gravitino is not thermalized in the ex- 
panding universe. Although it is out of the thermal bath, however, gravitinos are produced 
via the scattering processes of the thermal particles. In the expanding universe, evolution 
of the number density of the gravitino is governed by the following Boltzmann equation 

^^^^ + 3Hny2 = (atot^^rel)^rad) (F-2) 

where (ctot'^^rci) is the "thermally averaged" total cross section (times relative velocity), 
and nrad = ^^T^- With the Lagrangian ()F.1|) . the gravitino production cross sections 
are typically O {gf /327!-M^), where gi denotes the gauge coupling constant for SU{3)c, 
SU{2)l and U{1)y for z = 3, 2 and 1, respectively. 

Approximated formula for the primordial gravitino abundance is obtained by inte- 
grating Eq. ()F.2|) from the highest temperature in the radiation-dominated epoch, which 
approximately corresponds to the "reheating temperature" Tr after the infiation. Using 
the fact that the gravitino production cross section has very weak dependence on T, Eq. 
(IF.2|1 can be easily solved. In order to parameterize the primordial abundance of the 
gravitino, it is convenient to define the "yield variable" as ^^"^ 

Y,/. ^ ^, (F.3) 
s 

where s is the entropy density which is given by 

. = %9.sT\ (F.4) 

where g^^s is the effective number of the massless degrees of freedom. For the particle 
content of the minimal supersymmetric standard model (MSSM), g^^s = 228.75 for tem- 
perature much higher than the masses of the superparticles while g^,s = 43/11 for T -C mg. 
As we see just below, gravitino production processes are effective when T ~ Tr and hence, 
when T <^ Tr, becomes constant (as far as the total entropy in the comoving volume 
is conserved) . Then, integrating Eq. ()F.2j) from T = Tr to T ^ Tr, we obtain 



^ 45C(3) 



(o-tot'Wrol) ^rad 

g*sH 



(F.5) 



T=Tn 



'^^''The yield variable used in this paper differs from that in |5] which defined Y^^^^ = 713/2 /?^rad• We 
use ¥3/2 defined in Eq. ljF.3|) since, with this definition, Y3/2 is independent of time for the temperature 
T ^ Tr as far as there is no entropy production. Notice that, for T <C me, Y^J^^^ ~ 14Y3/2- 
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The important point is that the primordial abundance of the gravitino is (approximately) 
proportional to the reheating temperature after the inflation. Thus, gravitinos are over- 
produced in many cases if the reheating temperature is too high. 

For the detailed calculation of ^3/2, however, the definition of the "reheating tem- 
perature" used in the above calculation is quite ambiguous since, in realistic models of 
slow-roll inflation, there exists a period when the universe is dominated by the oscillating 
energy of the infiaton field before the radiation-dominated epoch. (We call this period 
as the inflaton-dominated period.) Due to the decay of the infiaton field, the inflaton- 
dominated epoch evolves into the radiation-dominated epoch. The transition from the 
inflaton-dominated epoch to the radiation-dominated epoch occurs when the expansion 
rate of the universe is comparable to the decay rate of the infiaton. In this paper, we 
define the "reheating temperature" as 

T^^(^^M^TfS\ (F.6) 

where (7* = (7*5. (In estimating Tr, we use = 228.75.) Notice that Tr given above 
is the same as the reheating temperature derived with the approximation such that X 
instantaneously decays when the relation 3H = Finf is realized. 

Then, in order to take account of the effects of the gravitino production in the inflaton- 
dominated epoch, we numerically solve the Boltzmann equations. The Boltzmann equation 
for the gravitino abundance is given in Eq. ()F.2|) . while evolutions of the energy densities 
of the radiation p^ad and infiaton p^^f are governed by the following Boltzmann equations 

^P''"^ + AHp,^ = FinfPinf, (F.7) 



dt 

dpini 
dt 



+ Sifpinf = -FinfPinf, (F.8) 



where Fjnf is the decay rate of the infiaton, and prad is related to the background temper- 
ature as 

Prad = ^9.T\ (F.9) 

The thermally averaged cross section for the gravitino production is calculated in 
Ref. [H2|- For SU{N) super Yang-Mills model with nf pairs of fundamental and 

'^^^In the gravitino production cross section, infrared singularities arise due to the t- and u-channel 
exchanges of the gauge bosons. In the previous calculation of (crtotWroi) given in such singularities 
were smeared by introducing a cut-off parameter for the scattering angle while, in |82| . it is treated by 
properly taking account of the thermal mass of the gauge bosons. Numerically, difference of the thermally 
averaged cross sections given in Refs. (9j and (82| is at 10 % level for most of the parameter space. 
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anti-fundamental chiral superfields, we obtain 



(O"totfrel) 



1 + 



X 



71 



327rM2 



C(3) 



{ [ln(TVmJ_th) + 0.3224] (A^ + m) + 0.5781nf} , (F.IO) 



where rrig is the gaugino mass and rug^t^ is the thermal mass of the gauge boson which is 
given as 



m 



g,th 



6 



g\N + n,)T'. 



(F.ii: 



As the gaugino mass becomes larger, more gravitinos are produced. In our calculation, 
we calculate the gravitino abundance in the limit of rrig ^ to derive a conservative con- 
straint. Thus, in our calculation, the gravitino abundance is underestimated by (9(10) % 
when the gaugino masses are comparable to the gravitino mass. If the gauginos are much 
heavier than the gravitino, upper bound on Tr is approximately obtained by rescaling 



the results by the factor 



m 
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/3m 



■3/2J 



using the fact that the gravitino production is 



dominated by processes related to the SU{3)c gauge fields. 

We follow the evolutions of 7x3/2, Prad, and pinf from the period with H ^ Finf to 
H <^ Finf by numerically solving the Boltzmann equations ()F.2|) . ()F.7|) and ()F.8|) . In 
calculating (crtot'^rei) , we sum over the contributions from all the MSSM gauge groups, 
SU{3)c, SU{2)l and U{1)y- Then, we calculate Y3/2 aX H <^ Fjnf (i.e., T ^ Tr). As we 
mentioned, F3/2 is approximately proportional to Tr; more precisely, we found that the 
resultant gravitino abundance is well approximated by the following formula 



3/2 



1.9 X 10"^2 
Tr 



X 



IQio GeV 



1 + 0.045 In 



T 



R 



1010 GeV 



1 - 0.028 In 



T 



R 



1010 GeV 



(F.12) 
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